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\  Abstract 
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This  thesis  describes  the  conversion  of  Program  EAGLE  -  Numerical 
Grid  Generation  System  to  ANSI  FORTRAN77  standards;  The  Development  of 
the  NAMELIST  input  emulator  for  Program  EAGLE;  and  the  modifications  and 
additions  allowing  interactive/non-interactive  execution  of  Program 
EAGLE.  All  the  work  being  done  on  an  IRIS  4D/70GT  computer  graphics 
workstation.  Description  of  additional  user  supplied  inputs  for 
improved  input  flexibility  is  also  given  as  well  as  a  discussion  of  the 
logic  implemented  in  the  routines  used  by  the  NAMELIST  input  emulator. 

Results  are  presented  for  three  generic  airframes  demonstrating  the 
capabilities  of  Program  EAGLE  as  an  IRIS  4D/70GT  computer  graphics  work¬ 
station.  Results  discussing  the  porting  of  the  standardized  version  of 
the  EAGLE  code  to  a  SUN  H/280  computer  system  are  also  presented. 
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Chapter  1 


Introduction 

With  the  advent  of  large-scale  computing  machines  (Cray  X-MP,  CRAY 
2,  etc.),  numerical  techniques  in  computational  fluid  dynamics  (CFD)  are 
now  able  to  solve  for  the  flow  field  surrounding  complex  airframe  con¬ 
figurations1  .  Numerical  grid  generation  has  been  cited  as  a  major  pac- 
ing  item  for  realistic  aircraft /missile  applications^  and  enables 
researchers  to  discretize  the  domain  about  arbitrarily-shaped  geo¬ 
metries. 

Considerable  progress  has  been  made  in  surface  (boundary)  and  grid 
(mesh)  genera tion^"^.  These  advances  allow  for  the  application  of 
multi-block  elliptic  grid  generation  techniques  to  complex  aerospace 
vehicles  as  well  as  to  other  complex  configurations  from  other  disci¬ 
plines  (i.e.  electromagnetics,  hydrodynamics,  etc.). 

In  1987,  the  CFD  community  saw  the  introduction  of  a  generalized 
three-dimensional  arbitrary  geometry  grid  generation  code  called  Program 
EAGLE  -  Numerical  Grid  Generation  System^ .  Program  EAGLE  is  a  compos¬ 
ite  (multi-block),  algebraic  and  elliptic  grid  generation  system  de¬ 
signed  to  discretize  the  domain  in  or  around  any  arbitrarily  shaped 
three-dimensional  region.  The  code  combines  a  three-dimensional,  sur¬ 
face  generation  scheme  with  a  multi-block,  three-dimensional 
boundary-conforming  elliptic  grid  generation  scheme. 
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The  surface  generation  system  of  Program  EAGLE  serves  as  a  front- 
end  to  the  grid  generation  system.  The  surface  generation  routine  de¬ 
velops  surfaces  (or  curves  in  2D)  to  be  input  to  the  grid  generation 
routine  as  segments  of  the  boundary  of  the  region  within  which  the  grid 
is  to  be  constructed^’^. 

The  grid  generation  system  is  a  general  two  or  three-dimensional 
algebraic  and  elliptic  grid  generation  routine  based  on  a  block  struc¬ 
ture,  which  allows  any  number  of  blocks  to  be  used  to  fill  an  arbitrary 
two  or  three-dimensional  region*0.  Any  block  can  be  linked  to  any  other 
block  (or  to  itself)  with  complete  (or  lesser)  continuity  across  the 
block  interfaces  as  specified  by  input.  This  routine  uses  an  elliptic 
generation  system  with  automatic  evaluation  of  control  functions.  These 
are  evaluated  either  directly  from  the  initial  algebraic  grid  and  then 
smoothed,  or  interpolated  from  the  boundary  point  distributions'^ *8»9,10. 

Other  features  of  the  surface  generation  routine  and  the  grid  gen¬ 
eration  routine  are  discussed  in  detail  in  References  9,  10,  11,  and  12. 

Program  EAGLE  was  initially  developed  to  execute  on  large-scale 
computing  machines  (supercomputers)  in  a  non-interactive  environment^’®. 
Since  its  introduction,  though,  interest  has  increased  in  developing 
interactive  grid  generation  programs  to  execute  on  engineering  computer 
workstations,  as  well  as  mini-super  computers  and  also  supercomputers 
(such  as  the  Cray  2).  The  increased  power  and  capabilities  of  these 
computer  systems  have  further  sparked  the  interest  of  computational 
fluid  dynamic  is ts  and  computational  aerodynamic  is ts  in  developing  nu¬ 
merical  grid  generation  techniques  in  an  interactive  environment. 
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Thus,  the  objective  of  this  work  is  to  develop  an  interactive/non- 
interactive,  portable,  user-oriented  version  of  the  EAGLE  code  on  a 
state-of-the-art  computer  workstation,  and  to  demonstrate  the  capabilities 
of  the  interactive  version  of  the  code  using  complex  aerospace  configu¬ 
rations.  A  design  criterion  was  that  changes  made  in  the  code  be  trans¬ 
parent  to  the  user  and  not  require  any  changes  of  the  inputs  to  the 
EAGLE  code. 

The  purpose  of  this  thesis  is  to  describe  this  effort  to 
knowledgeable  users  of  Program  EAGLE  interested  in  using  the 
standardized  interactive  version  of  the  EAGLE  code.  This  thesis  can 
and  should  be  considered  as  a  user's  manual  of  the  standardized 
interactive  version  of  the  EAGLE  code  and  as  a  supplement  to  the  user's 
manual  of  Program  EAGLE  in  Reference  9. 
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Chapter  2 


Developments 

Program  EAGLE  -  Numerical  Grid  Generation  System  was  originally  de¬ 
veloped  for  the  Cray  1  and  Cray  X-MP  supercomputing  systems'^.  The  most 
recent  version  of  Program  EAGLE  has  continued  to  be  developed  on  the 
Cray  X-MP  and  Cray  2  computer  systems^.  The  code  is  written  in  the  FOR¬ 
TRAN  language  available  on  these  systems  and  uses  the  NAMELIST  input  ex¬ 
tension  also  available  on  these  systems^’ ^ . 

Although  Program  EAGLE  is  written  in  the  FORTRAN  language,  it  does 
not  completely  follow  the  FORTRAN77  standards13*14.  Following  these 
standards  will  allow  codes  to  be  ported  over  to  other  computer  systems 
having  the  FORTRAN77  language  compiler  in  a  much  easier  fashion,  i. e. 
codes  can  compile  and  execute  with  little  or  no  modifications  on  differ¬ 
ent  computer  systems. 

The  other  problem  of  Program  EAGLE  is  that  it  uses  the  Cray  NAME- 
LIST  Input  extension1^.  This  feature  is  very  beneficial  in  that  it  al¬ 
lows  users  to  supply  the  inputs  into  the  code  by  specifying  the 
appropriate  variables  with  their  resnc'tive  values  in  any  order  in  the 
input  line.  This  makes  the  list  of  inputs  much  easier  to  understand  and 
reduces  the  reliance  on  the  user 's  manual  that  would  be  required  if  the 
inputs  were  structured  (i.e. ,  list  -  directed  and/or  formatted  inputs). 
As  mentioned  earlier,  NAMELIST  input  is  an  extension  and  not  a  standard 
of  F0RTRAN77.  Therefore,  not  all  computer  systems  carrying  the 
FORTRAN77  compiler  have  this  feature  or  extension  (for  example,  computer 
workstations,  such  as  the  IRIS-4D/70GT,  and  others). 
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To  be  able  to  satisfy  the  stated  objectives,  then,  requires  the 
EAGLE  code  to  be  converted  to  the  F0RTRAN77  language  standards.  This 
also  means  either  eliminating  the  NAMELIST  input  feature  or  writing  a 
NAMELIST  input  emulator  ( in  FORTRAN77  standard) .  To  ensure  the  changes 
made  in  Program  EAGLE  are  transparent  to  the  user  required  that  a 
NAMELIST  input  emulator  be  written. 

2.1  -  Conversion 

Converting  Program  EAGLE  to  F0RTRAN77  standards^  required  that 
both  routines  of  the  numerical  grid  generation  system  (i.e.,  surface  and 
grid  generation  routines)  be  converted  independently  and  compiled  cor¬ 
rectly.  Since  both  routines  have  similar  logic,  converting  the  first  one 
was  the  most  time  consuming  and  most  difficult,  the  reason  being  that 
knowledge  and  experience  had  to  be  gained  of  the  logic  and  the  routine's 
operation.  Once  this  knowledge  and  experience  was  gained,  the  second 
routine  was  converted  quite  easily  and  quickly  to  FORTRAN77  standards. 

The  most  obvious  first  step  was  to  eliminate  the  Cray  NAMELIST  in¬ 
put  extension  within  the  codes.  The  next  step  required  changes  in  using 
Hollerith  (or  character)  strings.  FORTRAN77  standards  do  not  accept 
Hollerith  strings  but  does  use  character  strings  instead.  This  was  the 
major  portion  of  the  conversion. 

Converting  from  Hollerith  to  character  manipulation  required  all 
character,  integer  and  real  variables  to  have  only  character,  integer 
and  real  values,  respectively.  This  was  accomplished  by  specifying 
variables  of  type  CHARACTER  to  be  character  variables,  variables  of  type 
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INTEGER  to  be  integer  variables  and  variables  of  type  REAL  to  be  real 
variables.  To  allow  the  user  the  flexibility  of  specifying  character 
values  for  integer  and  real  variables,  the  character  values  are  con¬ 
verted  to  unique  integer  (or  real,  respectively)  values.  These  values 
are  set  in  PARAMETER  or  DATA  statements  in  the  main  portion  of  each  rou¬ 
tine.  Note  that  these  unique  numbers  can  and  may  need  to  be  changed  to 
accommodate  the  computer  system  on  which  these  routines  are  to  execute. 

FORTRAN77  does  not  allow  character  values  to  be  embedded  in  common 
blocks  where  integer  and/or  real  variables  exist  within  the  same  common 
block.  Therefore,  all  character  variables  were  removed  from  these  com¬ 
mon  blocks  and  passed  to  all  the  necessary  routines  through  the  CALL  and 
SUBROUTINE  statements.  Separate  common  blocks  for  these  character  vari¬ 
ables  coUld  have  been  included.  But,  as  the  development  of  the  NAME- 
LIST  input  emulator  was  initiated,  it  was  easier  to  include  the  charac¬ 
ter  variables  in  the  CALL  and  SUBROUTINE  statements  instead  of  common 
blocks. 

With  the  major  part  of  the  conversion  completed,  other  smaller 
modifications  were  still  needed  to  allow  proper  compilation  and  execu¬ 
tion.  One  of  these  included  the  removal  of  various  DATA  statements  for 
variables  residing  in  common  blocks.  FORTRAN77  does  not  allow  initi¬ 
alization  of  variables  residing  in  common  blocks  by  DATA  statements  un¬ 
less  the  BLOCK  DATA  subprogram  feature  is  used.  Again,  because  of  the 
development  of  the  NAMELIST  input  emulator,  it  was  easier  to  initialize 
these  variables  using  executable  statements  at  the  beginning  of  the  pro¬ 
gram  instead  of  using  the  BLOCK  DATA  subprogram  feature  of  FORTRAN77. 
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Other  changes  to  the  routines  were  made  to  follow  suggestions  and 
guidelines  given  in  Reference  14.  One  of  these  is  to  use  the  generic 
names  of  intrinsic  functions  wherever  possible.  Where  intrinsic  func¬ 
tions  are  used  within  the  routines,  the  names  of  these  functions  were 
changed  to  their  generic  names  from  their  specific  names  wherever  possi¬ 
ble. 

Another  change  included  the  opening  and  closing  of  files  within  the 
routines  instead  of  allowing  the  system  to  open  and  close  them.  This 
gives  the  user  greater  flexibility  and  control  over  the  files  needed. 

After  all  these  changes  were  accomplished,  Program  EAGLE  compiled 
correctly  on  the  IRIS  4D/70GT  computer  graphics  workstation.  The  next 
steps  involved  the  development  of  the  interactive  logic  and  finally  the 
NAMELIST  input  emulator. 

2.2  Interactive  Logic 

The  logic  of  Program  EAGLE  is  such  that  converting  the  code  to  in¬ 
clude  interactive  processing  was  relatively  simple  with  minimal  changes. 

First,  a  status  indicator  (ISTAT)  was  included  to  indicate  the  er¬ 
ror  status  -  zero  (0)  for  no  error  (default,  of  course),  a  negative  one 
(-1)  for  a  "hard"  error  and  a  positive  one  (1)  for  a  "soft"  error.  A 
"hard"  error  occurs,  for  example,  when  the  problem  size  has  exceeded 
preset  dimensions,  thus  terminating  execution.  A  "soft"  error  occurs 
when  the  inputs  are  incorrect  or  missing.  The  routines  then  prompt  the 
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user  to  retype  the  input  line  correctly.  A  "soft"  error  behaves  just 
like  a  "hard"  error  in  the  non- interactive  or  batch  mode  (Reference  9 
has  details  on  the  errors  that  may  occur). 

This  variable  ISTAT  (either  set  to  1  or  -1 )  with  a  GO  TO  statement 
replaces  the  STOP  statements  found  in  the  original  version  after  an  er¬ 
ror  has  been  determined.  This  allows  execution  to  return  to  a  central 
location  within  the  routines  and  determines  the  next  course  of  action 
depending  on  the  value  of  ISTAT. 

Lastly,  the  variable  IBATCH  determines  whether  interactive  or  non¬ 
interactive  processing  is  to  take  place.  A  value  of  zero  (default)  in¬ 
dicates  interactive  execution,  while  a  value  of  one  (1)  indicates  non¬ 
interactive  execution  of  the  code.  Note  that  if  IBATCH  is  set  to  one 
(1)  and  interactive  execution  is  attempted,  the  program  will  behave  as 
if  it  is  executing  non-interactively.  Therefore,  any  errors  (hard  or 
soft)  will  terminate  execution.  IBATCH  also  controls  the  writing  of  the 
prompt  to  standard  output.  A  value  of  one  (1)  will  not  write  the 
prompt. 

2.3  Generalized  Parser /NAMELIST  Input  Emulator 

The  FORTRAN77  does  not  have  a  NAMELIST  input  feature  included  as  a 
standard.  NAMELIST  input  exists,  only  as  an  extension  available  on  some 
computer  systems.  Also,  NAMELIST  input  may  not  be  the  same  between  sys¬ 
tems  having  this  feature,  thereby,  making  the  original  EAGLE  code  dif¬ 
ficult  to  port  over  to  various  computer  systems,  especially  systems 
which  do  not  have  the  NAMELIST  input  feature  available.  To  be  able  to 
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port  Program  EAGLE  to  different  computer  systems,  and  enable  users  to 
use  their  previously  developed  input  lists,  a  generalized  parser  was 
written  to  emulate  the  Cray  NAMELIST  input  feature^. 

The  inclusion  of  the  emulator  into  the  EAGLE  code  (i.e,  into  the 
surface  and  grid  generation  routines)  did  not  change,  alter  or  modify 
the  structure  of  the  code;  however,  additional  arrays  were  required.  The 
addition  of  these  arrays  does  not  increase  the  memory  requirements  sig¬ 
nificantly,  since  the  additional  arrays  are  single-dimension  arrays  that 
are  two  orders  of  magnitude  smaller  than  the  arrays  already  present  in 
the  code. 

The  inclusion  of  the  emulator  does  require  additions  into  the  call¬ 
ing  routines.  Although  not  as  simple  as  adding  the  NAMELIST  input  fea¬ 
ture  into  a  code  as  described  in  Reference  15,  the  additions  are  not 
complex.  These  additions  are  non-executable  statements  such  as  PARAME¬ 
TER,  DIMENSION,  EQUIVALENCE  and  DATA  statements. 

Other  additions  include  the  two  CALL  statements  to  the  parsing  rou¬ 
tines  of  course.  The  first  CALL  statement  calls  the  subroutine  RDSTRG 
which  reads  in  a  line  as  a  character  string  with  a  maximum  of  132  char¬ 
acters.  The  character  string  is  then  passed  into  the  calling  routine  to 
be  passed  to  subroutine  PARSER,  which  is  the  actual  NAMELIST  input  emu¬ 
lator. 

once  control  returns  to  the  calling  routine,  the  inputs  are  passed 
into  the  necessary  variables  by  use  of  EQUIVALENCE  statements.  The  fact 
that  EQUIVALENCE  statements  were  needed  to  communicate  from  the  parser 
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to  the  calling  routines  required  that  some  of  the  labeled  common  blocks 
be  removed  or  modified  due  to  the  restrictions  imposed  by  using  common 
blocks  and  EQUIVALENCE  statements. 

All  the  additions  and/or  modifications  made  can  be  identified  in 
the  routines  by  searching  for  the  line(s)  containing  the  date  in  which 
the  additions  (and/or  modifications)  were  done  followed  by  the  charac¬ 
ters  "<gus  AM>".  The  end  of  the  sections  with  the  additions  and/or 
modifications  is  signified  by  a  line  beginning  with  a  "C"  followed  by  a 
string  of  asterisks. 
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Chapter  3 


User  Inputs 

All  inputs  are  given  in  the  NAMELIST  input  format  of  the  form 
E$INPUT  ITEM  *  'operation',  quantity  *  value,  .  ..$ 

The  "E"  in  column  1  indicates  to  print  (or  echo)  the  input  line  onto  the 
output  file.  If  the  "E"  is  omitted,  then  column  1  must  either  be  blank 
or  have  a  "C"  (meaning  to  ignore  this  line). 

The  first  dollar  sign  ($)  indicates  the  beginning  of  the  input 
line,  while  the  second  one  signifies  the  end  of  the  input  line.  The 
name  of  the  NAMELIST  follows  the  first  dollar  sign,  with  a  required 
space  after  the  name.  INPUT  is  the  name  of  the  NAMELIST  in  this  exam¬ 
ple.  The  names  used  for  the  NAMELIST  in  this  version  of  the  EAGLE  code 
are  SINPUT  -  for  inputs  to  the  surface  generation  routine  -  and  GINPUT 
and  GOUTPUT  -  for  inputs  to  the  grid  generation  routine. 

ITEM  ■  'operation'  designates  the  desired  operation,  while  the  val¬ 
ues  relevant  to  this  operation  follow  (i.e. ,  quantity  -  value).  In 
this  specification,  "quantity"  refers  to  the  variable  or  name  of  the 
input  quantity,  and  "value"  is  its  value. 

Values  for  arrays  can  be  given  in  one  of  two  ways.  They  are 
...,  quantity  =*  value,  value,  value,... 
or 

...»  quantity  =  N*value, . . . 

where  N  is  the  number  of  values.  Also,  a  value  for  a  particular  element 
of  an  array  can  be  specified  by  the  use  of  the  notation 

quantity  (M)  ■  value 
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where  M  indicates  the  element  of  the  array. 

Quantities  (variables)  of  type  CHARACTER  can  only  be  given  charac¬ 
ter  values  (with  a  maximum  of  64  characters).  Quantities  of  type  INTE¬ 
GER  and  REAL  can  be  given  integer  and  real  values,  respectively,  as  well 
as  specific  character  values  allowed  by  the  code  (see  Reference  9,  Vols 
II  and  III,  for  these  specific  values).  Presently,  the  NAMELIST  input 
emulator  is  not  able  to  accept  exponential  notation  (i.e.,  1.0E-04)  as 
input. 

For  the  purposes  of  clarity,  an  input  line  can  consist  of  any  num¬ 
ber  of  (continuation)  lines  where  each  line  has  a  maximum  of  132  charac¬ 
ters  (columns)  and  the  first  column  must  be  blank.  This  enables  an 
input  line  to  have  as  many  lines  as  needed,  as  long  as  the  first  line 
has  the  dollar  sign  in  the  second  column  and  the  last  line  has  the  sec¬ 
ond  dollar  sign  signifying  the  end  of  the  input  line.  Also,  all  inputs 
must  be  in  capital  letters  unless  specifying  a  filename.  The  filename 
in  quotes  can  be  lower  or  upper  case  letters.  And  finally,  single  quotes 
must  be  used  for  all  character  input  values. 

Adding  the  capability  of  interactive  execution  has  allowed  the  user 
to  have  greater  control  over  the  execution  of  the  code  by  the  addition 
of  three  new  operations. 

The  first  operation  allows  the  user  to  save  a  session.  The  input 
line  looks  like  the  following  example: 

E$SINPUT  ITEM  =  'SAVE',  FILNAM  =  'filename'  $ 
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where  "filename"  is  the  name  of  a  file  consisting  of  no  more  than  64 
characters  (includes  alpha-numeric  characters  and  symbols).  This  input 
line  should  be  the  first  line  of  input.  This  feature  only  operates  in 
the  interactive  mode.  Each  input  line  following  this  line  will  be  writ¬ 
ten  to  the  specified  file  (i.e. ,  "filename"). 

The  second  operation  allows  the  user  to  read  in  a  previously  saved 
session  or  a  previously  developed  input  file.  This  input  line  looks 
like  the  following  example: 

E$SINPUT  ITEM  *  *  READ ' ,  FILNAM  -  ’filename’^ 

Once  again  "filename"  is  the  name  of  the  file  to  be  read  in  by  the  code. 
This  input  line  will  cause  the  inputs  to  be  read  in  from  a  file  by  the 
name  of  "filename"  instead  of  reading  the  inputs  from  the  screen  (or 
standard  input) . 

The  third  operation  available  allows  the  user  to  terminate  execu¬ 
tion  when  desired  without  necessarily  providing  all  the  required  inputs 
for  a  particular  problem.  The  input  line  looks  like  the  following: 

E$SINPUT  ITEM  -  ’STOP* $ 

This  will  terminate  execution  and  save  all  the  necessary  files  if  this 
line  is  not  encountered  within  an  input  file.  In  other  words,  the  exis¬ 
tence  of  this  line  in  an  input  file  from  which  the  inputs  are  being  read 
will  terminate  reading  the  inputs  from  this  file,  and  then  will  expect 
the  rest  of  the  inputs  to  come  from  the  screen  (or  standard  input). 

Other  options  added  to  Program  EAGLE  give  the  user  greater  flexi¬ 
bility  in  choosing  the  type  of  output  desired.  Within  the  surface  and 
grid  generation  routines  the  user  has  the  option  of  specifying  the 
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filenames  on  input,  or  using  the  default  filenames  of  the  files  to  be 
written  or  read  in.  This  can  be  done  by  including  the  following  in  the 
input  line: 

FILNAM  -  •filename',... 

where  "filename"  is  the  name  of  the  file  to  be  written  or  read  in. 

For  example,  an  input  line  to  the  surface  generation  routine  may 
look  like  this: 

E$SINPUT  ITEM  -  'COMBINE',  COREIN  -  1,  -5, 

FILEOUT  -  1,  FILNAM  -  ' BLOCK 1 ' $ 

This  input  line  will  combine  the  segments  stored  in  cores  1  thru  5  and 
write  them  out  on  the  file  named  "BL0CK1".  The  grid  generation  routine 
can  then  read  in  this  file  with  the  following  example  input  line: 

E$GINPUT  ITEM  -  'FILE*,  FILE  -11,  FILNAM  -  ' BL0CK1 ',...$ 

Another  additional  feature  allows  the  user  to  select  the  output 
from  the  grid  generation  routine  in  the  form  needed  to  be  read  in  by  the 
NASA-AMES  graphics  program  called  PL0T3D.  To  select  this  feature,  the 
user  must  include  the  following  in  an  input  line: 

...,  OUTER  -  'PL0T3D',  ... 

This  option  is  generally  selected  when  the  file  on  which  the  grid  is  to 

be  written  is  specified.  For  example: 

E$GINPUT  ITEM  -  'STORE',  FILE  -  72,  FILNAM  =  'FILE72. fmt ' , 

OUTER  -  ' PL0T3D*  $ 

Selecting  the  output  feature  "PL0T3D"  will  write  formatted  records  to 
the  FORTRAN  file  72,  which  has  the  filename  »FILE72.fmt"  in  the  form 
needed  by  PLQT3D. 
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OUTER  *  'PLOT3D'  allows  two  options.  One  option  allows  the  user  to 
write  to  the  file  without  a  blanking  array  (default),  and  the  other  op* 
tion  will  allow  the  user  to  write  to  the  file  with  a  blanking  array.  To 
select  the  second  option  the  user  must  specify  the  operation  ITEM  - 
•BLANK'  and  its  associated  quantities  and  values  in  the  output  phase  of 
the  grid  generation  routine.  For  example,  an  input  file  to  the  grid 
generation  routine  selecting  the  PL0T3D  option  with  blanking  may  look 
like  the  following: 


E$GINPUT  ITEM  -  'STORE' ,  FILE  *  72,  FILNAM  -  ' FILE72. fmt ' , 

OUTER  -  'PL0T3D'  $ 


E$GINPUT  ITEM  -  'END'  $ 

E$GOUTPUT  ITEM  -  'BLANK',  BLOCK  -  1,  START  =■  1,1,1, 
END  -  65,1 ,31 .  VALUE  -  2  $ 


E$G0UTPUT  ITEM  =  'END'  $ 

In  this  input  list,  file  72  will  have  the  filename  of  "FILE72. fmt" ,  and 
the  records  will  be  written  formatted  in  the  form  needed  to  be  read  in 
by  the  program  PL0T3D.  With  the  operation  ITEM  ■  'BLANK'  selected, 
file  72  will  also  contain  an  integer  array  called  IBLANK  specifying 
which  points  are  to  be  'blanked'.  In  this  example,  in  block  1  the 
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points  starting  at  1,1,1  to  65,1,31  will  have  an  integer  value  of  2 
(VALUE-2)  stored  in  the  IBLANK  array.  Any  integer  value  besides  zero 
(0)  in  the  IBLANK  array  will  tell  the  PL0T3D  program  to  draw  that  point 
on  the  screen.  Any  integer  value  can  be  specified  in  the  quantity 
VALUE.  The  default  though  is  one  (1). 
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Chapter  H 


Generalized  Parser /NAM EL 1ST  Input  Emulator  Operation 

As  mentioned  earlier,  Program  EAGLE  was  originally  written  to  use 
the  Cray  NAMELIST  inout  extension^ **,  which  is  not  available  on  all 
computer  systems.  Therefore,  to  ensure  uniformity  of  inputs  on  all  com¬ 
puter  systems,  and  to  allow  input  lists  previously  developed  to  be  used 
without  making  changes  to  these  lists,  a  generalized  parser/NAMELIST 
input  emulator  was  written.  This  emulator  generally  follows  the  basic 
rules  of  use  of  the  NAMELIST  input  extension  specified  in  Reference  15. 
Several  basic  guidelines  will  be  stated  here,  however,  for  the  purpose 
of  clarity. 

All  upper-case  letters  must  be  used  except  when  specifying  file¬ 
names.  The  filenames  can  be  given  in  lower  or  upper-case  letters.  Sin¬ 
gle  quotes  must  be  used  everywhere. 

An  input  line  can  consist  of  any  number  of  lines,  where  the  maxi¬ 
mum  number  of  characters  per  line  is  132.  An  input  line  begins  with  a 
dollar  sign  ($)  in  column  two  (2),  and  is  terminated  with  a  second  dol¬ 
lar  sign  at  the  end  of  the  input  line,  regardless  of  how  many  lines  make 
up  the  input  line.  Each  line  after  the  first  line  of  the  input  line 
must  start  on  or  after  column  two  (2). 

Column  three  (3)  of  the  first  line  must  begin  with  the  name  of  the 
input  list  (the  NAMELIST  name),  and  a  blank  must  separate  the  name  with 
the  actual  beginning  of  the  inputs.  Blanks  can  be  used  throughout  to 
make  the  inputs  more  legible,  but  blanks  cannot  appear  within  variable 
names  or  between  single  quotes  when  specifying  character  values. 
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Column  one  (1)  must  be  blank  or  have  either  an  "E"  or  a  "C".  An 
"E"  in  column  one  (1)  indicates  to  write  the  input  line  to  standard  out¬ 
put.  A  "C"  in  column  one  (1)  means  to  ignore  this  line;  however,  this 
line  will  be  written  to  standard  output. 

^.1  -  Parameters,  Variables  and  Arrays 

The  operation  of  the  parser/NAMELIST  input  emulator  is  basically 
divided  into  two  major  routines.  The  first,  subroutine  RDSTRG,  reads  in 
a  character  string  while  the  second,  subroutine  PARSER,  takes  this 
string,  checks  it  and  then  parses  it  to  the  appropriate  variables. 

In  the  routine  RDSTRG  the  variable  STRING  is  the  character  string 
to  be  read  in,  with  the  actual  maximum  character  length  being  specified 
by  the  variable  ISTRNG.  The  variable  NSTRNG  provides  the  maximum  al¬ 
lowed  character  length  for  the  variable  STRING.  Note  that  NSTRNG  should 
be  less  than  or  equal  to  ISTRNG.  The  variable  IREAD  determines  whether 
to  read  from  a  file  (IREAD  *  1 )  or  from  standard  input  (default).  The 
variable  ISTAT  specifies  the  error  status  -  zero  (0)  means  no  error 
while  any  other  number  (say,  1  and  -1)  indicates  an  error  has  occurred. 
The  variable  LCT  is  the  actual  length  of  the  character  string  read  in. 
The  variables  STRING,  LCT  and  ISTAT  are  passed  back  to  the  calling  rou¬ 
tine. 

The  character  string  STRING  is  then  passed  from  the  calling  routine 
into  subroutine  PARSER  through  the  CALL  and  SUBROUTINE  statements  for 
parsing.  The  parameters  NVAR,  NCHAR,  NINTG  and  NREALS  are  used  to  di¬ 
mension  the  arrays  in  subroutine  PARSER.  NVAR  is  the  total  number  of 
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variables  the  user  can  supply  as  input.  This  parameter  is  the  sum  of 
the  other  three  parameters,  NCHAR,  NINTG  and  NREALS.  These  three  are 
the  total  number  of  characters,  integers  and  real  variables,  respect¬ 
ively,  that  the  user  can  supply  as  input. 

The  three  parameter i>  NCVAL,  NIVAL  and  NRVAL  are  the  total  number  of 
characters,  integer  and  real  values,  respectively,  that  can  be  supplied 
during  input.  The  difference  between  these  three  parameters  and  the 
previous  set  of  parameters  is  that  these  three  parameters  contain  the 
sum  total  of  all  the  elements  of  the  arrays  of  the  input  variables. 
NCVAL,  NIVAL  and  NRVAL  can  be  the  same  as  NCHAR,  NINTG  and  NREALS  if  all 
the  variables  which  can  be  supplied  during  input  contain  only  one  ele¬ 
ment  (i.e.,  all  the  variables  are  not  arrays). 

The  last  three  parameters  are  NDAT,  NIDAT,  and  NRDAT.  The  sum  of 
NIDAT  and  NRDAT  gives  NDAT  -  the  total  number  of  integer  (NIDAT)  and 
real  (NRDAT)  data  constants  set  in  the  calling  routine.  These  are  used 
for  the  purpose  of  specifying  character  values  for  integer  and  real 
variables.  The  parser  then  knows  to  convert  supplied  character  strings 
into  unique  integer  and  real  numbers. 

The  arrays  needed  by  subroutine  PARSER  are  NAMLIS,  NDIM1 ,  NDIM2, 
ICHR,  INTGR,  REALS,  CDATA,  IDATA,  RDAT,  and  CPOS,  IPOS,  RPOS.  The  char¬ 
acter  array  NAMLIS  is  dimensioned  to  NVAR  and  contains  the  list  of  input 
variable  names  organized  with  input  variables  of  type  CHARACTER  first, 
type  INTEGER  next,  and  type  REAL  last.  This  array  is  used  to  check  the 
input  and  see  that  the  supplied  variable  exists  in  the  NAMELIST. 
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The  integer  arrays  NDIM1  and  NDIM2  are  dimensioned  to  NVAR  and  hold 
the  lower  and  upper  bound  dimensions,  respectively,  of  all  the 
variables.  The  values  within  these  arrays  are  organized  in  the  same 
manner  discussed  above.  A  two-dimensional  array  will  have  the  upper 
bound  of  its  dimension  set  to  the  total  number  of  elements  within  that 
array.  These  arrays  are  used  to  ensure  that  no  more  than  the  required 
number  of  values  are  passed  into  their  respective  variables. 

The  arrays  ICHR,  INTGR  and  REALS  are  character,  integer  and  real 
arrays,  respectively,  which  hold  the  values  of  all  the  variables.  These 
arrays  are  then  passed  back  into  the  calling  routine  and  equivalenced  to 
the  appropriate  variables.  These  arrays  are  dimensioned  to  NCHAR, 
NINTG,  and  NREALS,  respectively. 

The  character  array  CDATA  contains  the  names  of  ths  constants  set 
in  the  calling  routine.  This  array  is  used  to  check  the  input  and  see 
that  the  supplied  character  inputs  to  integer  and  real  variables  are 
present  within  this  array.  The  array  CDATA  is  dimensioned  to  NDAT. 

The  arrays  IDAT  and  RDAT  are  integer  and  real  arrays,  dimensioned 
to  NIDAT  and  NRDAT,  respectively.  These  arrays  contain  the  unique  inte¬ 
ger  and  real  constants  which,  will  be  passed  to  the  appropriate  variables 
when  character  values  are  supplied  to  integer  and  real  variables. 

The  last  three  arrays  are  CPOS,  IPOS,  and  RPOS.  These  three  are 
integer  arrays  dimensioned  to  NCHAR,  NINTG  and  NREALS,  respectively.  The 
purpose  of  these  three  arrays  is  to  hold  actual  array  positions  within 
the  three  arrays  ICHR,  INTGR  and  REALS,  respectively.  In  other  words 
CPOS,  IPOS  and  RPOS  act  as  pointers  for  the  parser  to  place  the  values 
in  the  proper  position  in  the  respective  arrays. 
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The  last  set  of  arguments  needed  by  subroutine  PARSER  consists  of 
two  character  strings  -  STRING  and  NAMLST  -  and  several  integer  vari¬ 
ables.  The  first  character  string  called  STRING  is  passed  from  the  rou¬ 
tine  RDSTRG  to  the  calling  program  and  then  to  subroutine  PARSER.  The 
second  character  string,  NAMLST,  holds  the  name  of  the  NAMELIST  input. 
This  variable  can  be  an  array  holding  different  names  of  the  NAMELIST 
input.  This  is  demonstrated  in  the  grid  generation  routine  of  the  EAGLE 
code. 

The  integer  variable  NDS,  indicates  where  the  second  dollar  sign  is 
located  in  the  character  string  STRING.  If  none  is  found  then  NDS  is 
set  to  zero  (0),  the  default. 

The  variable  L  maintains  a  sum  total  of  the  lines  read  in  per  input 
line.  The  default  value  for  L  is  one  (1). 

The  integer  variable  INL  stores  the  number  of  the  last  input  vari¬ 
able  found  in  STRING.  INL  also  has  a  default  value  of  one  (1). 

The  next  variable,  ITEXT,  keeps  track  of  the  number  of  values  sup¬ 
plied  for  the  last  input  variable  in  STRING. 

The  variable  IECH0  is  a  print  (or  echo)  indicator  which  is  set  to 
one  (1)  when  an  "E"  is  present  in  column  1  of  the  first  line  of  the  in¬ 
put  line.  The  default  value  for  IECH0  is  zero  (0). 

The  variable  LS  and  LT  store  the  character  lengths  of  the  character 
strings  STRING  and  TEXT,  respectively.  (TEXT  stores  the  name  of  the 
input  variable.)  Both  have  default  values  of  one  (1). 

The  variable  NL  serves  as  a  search  pointer  indicating  where  in  the 
character  string  STRING  the  search  last  ended.  Its  default  value  is  one 
(1). 
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The  next  variable  NJ  stores  the  total  number  of  elements  within  an 
array.  This  variable  is  used  to  ensure  that  the  value  stored  in  ITEXT 
is  less  than  or  equal  to  NJ.  Again,  the  default  value  for  NJ  is  one 
(1). 

IEQUAL  indicates  whether  an  equal  sign  is  present.  A  value  of  one 
(1)  means  an  equal  sign  was  found,  while  a  value  of  zero  (0)  means  no 
equal  sign  is  present. 

The  penultimate  variable,  NUM,  indicates  to  multiply  the  given 
value  by  an  integer  number  when  the  value  for  NUM  is  one  (1). 

The  last  variable  in  the  argument  list  (ISTAT)  is  the  status  indi¬ 
cator  discussed  earlier. 

All  these  parameters,  arrays  and  variables  are  needed  by  the  gen¬ 
eralized  parser/NAMELIST  input  emulator  to  successfully  parse  the  in¬ 
puts.  Only  the  character  string  STRING  and  the  three  arrays  storing  all 
the  inputs  -  ICHR,  INTGR  and  REALS  arrays  -  are  used  by  the  calling  rou¬ 
tine. 

k.2  Program  Configuration/Set-Up 

To  be  able  to  use  this  emulator,  several  non-executable  statements 
must  be  added  to  the  calling  routine^  as  well  as  two  CALL  statements  to 
the  routines  which  read  in  the  lines  of  input  (subroutine  RDSTRG)  and 
which  parse  these  inputs  to  the  appropriate  variables  (subroutine 
PARSER) .  Appendix  A  has  a  sample  listing  of  all  the  statements  which 
need  to  be  added  to  the  calling  routine.  The  following  paragraphs  will 
discuss  the  needed  statements. 
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To  begin  with,  several  arrays  and  the  variable  STRING  must  be 
specified  as  type  CHARACTER.  The  character  length  of  STRING  must  match 
the  value  specified  for  the  variable  ISTRING  in  the  PARAMETER  statement. 
The  value  for  NSTRING  can  be  set  to  ISTRNG  or  to  a  lesser  value  as  evi¬ 
denced  in  Appendix  A. 

Next,  PARAMETER  statements  containing  the  values  for  the  variables 
NCHAR,  NINTG,  NREALS,  and  NCVAL,  NIVAL,  NRVAL  and  NRDAT  are  required. 
Also,  PARAMETER  statements  will  be  needed  to  initialize  the  variables 
NVAR  and  NDAT.  NVAR  is  defined  by  the  equation 

NVAR  -  NCHAR  +  NINTG  +  NREALS 
and  NDAT  is  defined  by  the  equation 

NDAT  -  NIDAT  +  NRDAT 

These  statements  set  up  the  dimensions  needed  by  the  arrays  used  by  the 
routine  PARSER, 

The  next  item  is  to  dimension  the  arrays  needed  by  the  parser.  The 
arrays  NAMLIS,  NDIM1  and  NDIM2  are  dimensioned  to  NVAR.  The  arrays 
ICHR,  INTGR  and  REALS  are  dimensioned  to  NCVAL,  NIVAL  and  NRVAL,  respec¬ 
tively.  The  arrays  CPOS,  IPOS  and  RPOS  are  dimensioned  to  NCHAR,  NINTG 
and  NREALS,  respectively.  The  last  three  arrays  CDATA,  IDAT  and  RDAT 
are  dimensioned  to  NDAT,  NIDAT  and  NRDAT,  respectively. 

After  the  DIMENSION  statements  comes  the  setting  up  of  the  communi¬ 
cation  between  the  three  arrays  returned  by  the  parser  -  ICHR,  INTGR  and 
REALS  -  and  the  variables  used  by  the  calling  routine.  This  is  accom¬ 
plished  by  using  EQUIVALENCE  statements.  One  set  of  these  statements 
passes  the  values  in  the  character  array  ICHR  into  the  appropriate  charac- 
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ter  variables.  The  next  set  passes  the  integer  values  in  the  ar¬ 
ray  INTGR  into  the  appropriate  integer  variables.  And  finally,  the 
third  set  of  statements  passes  the  real  values  in  the  array  REALS  into 
the  appropriate  real  variables. 

The  last  set  of  EQUIVALENCE  statements  passes  unique  integer  and 
real  numbers  to  the  parser  through  the  arrays  IDAT  and  RDAT,  respect¬ 
ively.  This  last  set  of  statements  is  required  only  if  character  val¬ 
ues  are  to  be  supplied  during  input  to  integer  and  real  variables. 
These  unique  numbers  should  be  set  either  in  PARAMETER  statements  or 
DATA  statements. 

The  last  set  of  statements  -  DATA  statements  -  initialize  the  ar¬ 
rays  NAMLIS,  NDIM1  ,  NDIM2,  CPOS,  IPOS  and  RPOS,  as  well  as  initializing 
the  unique  numbers  and  the  default  values  for  the  variables  needed  in 
subroutine  PARSER. 

After  completing  the  set-up  of  the  non-executable  statements,  the 
section  which  reads  the  inputs  and  then  parses  the  inputs  must  be  added. 
The  first  CALL  statement  is  to  the  routine  RDSTRG  which  reads  the  input 
(as  a  character  string),  checks  its  length  and  returns  the  character 
string  STRING.  The  second  CALL  statement  is  to  the  routine  PARSER  where 
the  character  string  STRING  is  passed  into  it,  and  the  input  values  are 
returned  in  the  three  arrays  ICHR,  INTGR  and  REALS.  The  character 
string  STRING  can  be  used  by  the  calling  routine  since  it  is  not  changed 
in  any  way.  These  CALL  statements  should  be  placed  after  initializing 
the  defaults  and  prior  to  testing  for  correct  input  values. 
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4,3  Subroutines  and  Functions 


Subroutine  PARSER  not  only  governs  the  parsing  of  the  inputs  but 
also  has  the  logic  necessary  to  emulate  the  Cray  NAMELIST  input  fea¬ 
ture1  This  logic  is  the  first  thing  which  the  character  string  STRING 
goes  through. 

Initially  the  search  pointer  NL  is  set  b;  using  the  name  of  the 
NAMELIST  (NAMLST),  and  the  character  length  of  STRING  is  also  determined 
and  placed  in  the  variable  LC. 

The  first  section  of  the  logic  ensures  that  column  one  (1)  is  blank 
or  has  either  an  "E"  (writes  the  character  string  STRING  to  standard 
output)  or  a  ”C"  (the  line  is  ignored  and  written  out  to  standard  out¬ 
put  and  execution  is  returned  to  the  calling  routine).  After  this  test, 
a  dollar  sign  ($)  is  searched  for  in  column  two  (2).  If  found  then  the 
search  continues  for  the  second  dollar  sign  indicating  the  end  of  the 
input  line.  If  this  second  dollar  sign  is  not  found  in  the  first  line 
(L-1)  of  the  input  line,  the  variable  LS  is  set  to  the  character  length 
of  STRING  (LC).  Other  subsequent  lines  will  be  tested  until  the  second 
dollar  sign  is  found,  thus  terminating  the  input  line. 

The  next  test  in  the  logic  looks  for  the  name  of  the  NAMELIST 
(NAMLST)  in  STRING  starting  from  column  three  (3)  up  to  the  column 
specified  by  the  search  pointer  NL.  After  satisfying  this  condition,  a 
blank  in  column  NL+1  is  searched  for,  and  then  NL  is  reset  to  its  proper 
position  in  STRING. 
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The  next  part  in  subroutine  PARSER  will  put  together  a  string  of 
characters  in  the  character  string  TEXT  and  then  check  it  against  the 
list  (NAMLIS)  to  determine  if  the  supplied  input  variable  is  a  valid 
input.  If  so,  the  variable  INL  specifies  which  input  variable,  and  sets 
the  maximum  number  of  elements  (NJ)  and  the  element  position  within  one 
of  the  three  arrays  (ICHR,  INTGR  or  REALS)  for  that  particular  input 
variable.  Finally  STRING,  LS,  NJ  and  one  of  the  three  arrays  are  passed 
into  one  of  the  three  routines  (FNDCHR,  FNDINT  or  FNDRL)  for  parsing, 
depending  on  the  value  in  INL. 

4.3.1  FNDCHR 

Subroutine  RNDCHR  searches  for  a  character  value  within  the  charac¬ 
ter  string  STRING,  starting  from  the  last  position  of  the  pointer  (L)  to 
the  end  of  the  string  (LS)  or  until  the  character  value  has  been  found, 
whichever  comes  first. 

The  routine  begins  by  searching  for  the  equal  sign  after  the  input 
variable.  Generally,  the  variable  L  will  point  to  a  position  in  STRING 
which  is  a  blank  or  has  an  equal  sign  (-)  initially.  If  the  equal  sign 
is  not  found  prior  to  reaching  the  end  of  the  string,  execution  will 
return  to  the  routine  calling  the  parser. 

After  finding  the  equal  sign,  the  character  value  is  found  in 
STRING  by  searching  for  a  set  of  single  quotes  (•).  Once  the  first  sin¬ 
gle  quote  is  found,  routine  CHRVAL  is  called  to  get  the  position  of  the 
second  single  quote  within  STRING,  thus  indicating  a  character  value  is 
present.  The  character  value  is  then  stored  in  the  array  IC,  providing 
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that  the  number  of  values  found  (ITEXT)  does  not  exceed  the  number  al¬ 
lowed  (NJ).  The  search  pointer  (L)  is  then  updated,  and  control  is  re¬ 
turned  to  the  beginning  of  the  routine  to  determine  if  more  values  are 
present. 

If,  after  finding  the  equal  sign,  a  number  followed  by  an  asterisk 
(indicating  multiplication)  is  found,  then  the  upcoming  character  value 
will  be  repeated  in  the  array  IC  the  number  of  times  specified  by  this 
number. 

This  number  must  be  an  integer  and  is  stored  as  a  character  string 
in  TEXT  to  be  later  converted  to  an  integer  number.  Then  the  search 
continues  for  the  asterisk  (*).  Finding  the  asterisk,  routine  CONVER  is 
called  to  convert  the  character  string  in  TEXT  to  a  number  which  is 
stored  as  an  integer  number  in  ITEXT.  If  ITEXT  is  less  than  or  equal  to 
the  value  in  NJ  then  the  variable  NUM  is  set  to  one  (1),  and  the  search 
continues  for  the  character  value  as  described  above.  This  time, 
though,  the  character  value  is  stored  ITEXT  times  in  the  array  IC. 

The  search  is  terminated  and  control  is  returned  to  the  parser  once 
a  comma  and  a  character  after  the  comma  is  found  not  satisfying  any  of 
the  conditions,  and  the  value  in  ITEXT  is  less  than  or  equal  to  the 
value  in  NJ. 

4.3.2  FNDINT,  FNDRL 

Subroutines  FNDINT  and  FNDRL  search  for  an  integer  or  real  value, 
respectively,  within  the  character  string  STRING,  starting  from  the  last 
position  (L)  to  the  end  of  the  string  (LS),  or  until  all  the  integer  or 
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real  values  have  been  found  for  the  specified  input  variable.  Both  rou¬ 
tines  function  similarly,  with  the  exception  that  subroutine  FNDRL  re¬ 
quires  searching  to  include  a  decimal  point  (.).  Also,  both  routines 
use  the  same  logic  used  in  subroutine  FNDCHR  to  determine  when  the  input 
values  are  being  specified.  In  other  words,  an  equal  sign  is  found 
first,  and  then  the  input  values  come  after  the  equal  sign. 

After  finding  the  equal  sign,  an  integer  (or  real)  value  is  found 
in  STRING  by  searching  for  a  character  between  the  number  zero  (0)  and 
nine  (9),  or  a  plus  (  +  )  or  minus  (-)  sign.  In  subroutine  FNDRL  the 
search  also  includes  a  period  (.).  Once  any  of  these  symbols  are  de¬ 
tected,  they  are  placed  in  the  character  string  TEXT  as  characters  until 
a  comma  (,)  is  found.  Then  TEXT  is  converted  to  an  integer  (or  real)  value 
by  the  routine  CONVER  and  stored  in  the  array  IC  if  the  number  of  input 
values  (ITEXT)  does  not  exceed  the  maximum  specified  by  the  variable  NJ. 

Another  piece  of  logic  was  inserted  to  allow  input  of  character 
values  to  integer  and  real  variables.  The  character  values  are  found  by 
searching  for  a  single  quote,  calling  subroutine  CHRVAL  to  determine 
the  position  of  the  second  single  quote,  and  then  converting  this  charac¬ 
ter  value  to  a  unique  integer  (or  real)  value  using  subroutine  CCI  (or 
CCR),  respectively.  The  unique  integer  and  real  values  must  be  preset 
in  the  routine  calling  subroutine  PARSER,  as  well  as  setting  up  the  ar¬ 
rays  needed  by  the  two  routines,  CCI  and  CCR  (see  Section  H.2).  After 
conversion  to  unique  integer  (or  real)  values,  they  are  then  stored  in 
the  array  IC  as  described  earlier. 
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As  discussed  in  Section  4.3.1,  if  an  asterisk  is  present  after  a 
number  has  been  specified,  then  the  upcoming  input  value  will  be  re¬ 
peated  in  the  array  IC  the  number  of  times  specified  by  the  given  num¬ 
ber. 

The  search  is  terminated  and  control  is  returned  to  the  parser  as 
described  in  Section  4.3.1. 

4.3.3  CHRVAL 

Subroutine  CHRVAL  finds  a  character  value  in  the  character  string 
STRING  by  returning  the  position  of  the  second  single  quote  in  STRING. 

The  first  single  quote  is  found  prior  to  calling  this  routine,  and 
the  position  is  passed  into  CHRVAL  through  the  variable  I,  The  second 
single  quote  is  searched  for  in  STRING,  and  its  position  is  stored  in 
the  variable  IQ2.  If  no  second  quote  is  found,  IQ2  will  have  a  value  of 
zero  (0).  The  character  value  is  then  checked  for  embedded  blanks 
between  the  single  quotes,  and,  finally,  the  position  of  the  second  sin¬ 
gle  quote  is  returned  to  the  calling  routine  in  the  variable  IQ2, 


4.3.4  CONVER 


Subroutine  CONVER  converts  a  character  representation  of  an  integer 
or  real  number  to  an  actual  real  number  (excluding  exponential  represen¬ 
tation  of  a  real  number). 
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The  character  representation  comes  from  the  calling  routine  in  the 
character  string  TEXT.  TEXT  is  first  checked  to  determine  if  a  minus 
(-)  or  a  plus  (+)  sign  is  present.  If  a  minus  sign  is  found  then  the 
variable  SIG  is  given  the  value  of  -1,  while  if  a  plus  sign  is  found  SIG 
will  have  a  value  of  1.  The  beginning  (NSTR)  of  the  string  TEXT  is  then 
set  to  one  plus  the  present  point.  If  a  decimal  point  is  found,  the 
position  of  the  decimal  point  is  returned  in  the  variable  NDEC. 

Each  character  in  TEXT  is  converted  to  a  real  number  by  subtracting 
the  assigned  value  of  the  ANSI  representation  of  zero  (0)  from  the  as¬ 
signed  value  of  the  ANSI  representation  of  each  character  in  TEXT.  All 
the  numbers  are  then  added  to  obtain  the  value  stored  in  the  variable 
VAR.  The  result  in  VAR  is  a  real  number  which  can  always  be  converted 
to  an  integer  when  necessary  in  the  calling  routine. 

4.3.5  CCI,  CCR 

Subroutines  CCI  and  CCR  convert  character  values  passed  in  the 
character  string  TEXT  into  unique  integer  and  real  values,  respectively. 
These  unique  values  must  have  been  specified  in  the  routine  calling  sub¬ 
routine  PARSER  and  passed  through  the  arrays  IDAT  and  RDAT.  The  array 
CDATA  contains  the  list  of  character  values  to  be  specified  during  input 
for  integer  and  real  variables. 
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TEXT  is  first  checked  against  the  list  in  CDATA  to  ensure  that  the 
character  value  within  TEXT  is  allowed.  Then  the  position  of  the  char¬ 
acter  value  within  the  array  CDATA  is  checked  to  ensure  that  the  correct 
unique  number  is  stored  in  the  variable  IV AR  (VAR  in  subroutine  CCR) . 
IVAR  is  then  passed  to  the  calling  routine. 

M.3.6  ELEMNT 

Subroutine  ELEMNT  converts  a  character  representation  of  a  number 
within  a  set  of  parenthesis  to  an  integer  number  indicating  the  position 
in  the  array  specified  during  input. 

The  character  string  TEXT  contains  the  input  variable  (along  with 
the  parentheses  and  number)  in  character  representation.  The  positions 
of  the  left  and  right  parentheses  (NL  and  NR,  respectively)  within  TEXT 
are  determined.  Then,  between  the  parentheses,  the  character  represen¬ 
tation  for  the  number  is  stored  in  the  character  string  NUM.  NUM  is 
then  passed  into  the  routine  CONVER  for  conversion  to  an  integer  value 
indicating  which  element  in  the  array.  The  value  is  then  returned  in 
the  variable  IVAR. 

4.3.7  ILEN,  ISIZ 

Function  ILEN  determines  the  actual  number  of  characters  in  the 
character  string  TEXT. 

Function  ISIZ  counts  between  the  two  numbers  N1 ,  N2.  The  sum  total 
is  returned  in  ISIZ. 
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Chapter  5 


Results 


5.1  Applications 

Three  generic  configurations  were  used  to  demonstrate  Program  EA¬ 
GLE'S  capabilities  on  an  IRIS  40/70GT  computer  graphics  workstation.  The 
first  of  the  configurations  is  an  elliptic  cross-section  missile^*1^1^. 
The  second  is  an  ogive-cylinder-ogive  with  fins^'1®.  And  the  third  is 
an  elliptic  airframe  with  horizontal  and  vertical  control  surfaces®'1^. 

Each  configuration  was  generated  interactively  with  no  other  users 
on  the  system.  Also,  each  configuration  generated  used  the  option  in 
the  EAGLE  code  to  execute  as  .if  it  were  executing  on  a  Cray  X-MP  com¬ 
puter  system  with  solid  state  desk  device  (SSD).  In  other. words,  files 
which  are  assigned  to  SSD  on  a  Cray  X-MP  with  SSD  are  now  used  as  regu¬ 
lar  files.  Therefore,  CPU  times  stated  may  consist  of  up  to  80?  I/O 
time  and  not  actual  computation.  This  statement  is  made  based  on  simi¬ 
lar  runs  using  SSD  on  the  Cray  X-MP  computer  system1 

5.1.1  Elliptic  Missile 

This  configuration  has  an  elliptic  -  3:1  -  cross-section  throughout 
the  body  described  by  a  parabolic  curve  from  the  nose  to  the  aft  end  of 
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the  body.  Extending  from  the  end  is  an  attached  sting.  Figure  1  shows 
the  body  of  this  configuration. 

A  C-type  system  is  used  to  discretiize  the  domain  on  and  around  the 
body.  This  field  grid  is  described  by  two  (2)  blocks  (or  grids),  each 
block  containing  one-half  of  the  configuration.  Each  block  has  the  di¬ 
mensions  of  65x18x31  points  -  65  points  in  the  axial  direction,  18 
points  in  the  radial  direction  and  31  points  in  the  circumferential  di¬ 
rection  -  for  a  total  of  725*10  points  (36270  points  per  block). 

Figures  2a  through  3b  display  side  and  front  views  of  the  algebraic 
grid.  Figures  4a  and  4b  show  a  perspective  view  of  the  entire  grid  and 
a  closeup  view  of  the  grid  near  the  body.  The  surface  and  field  grids 
shown  in  these  figures  are  the  same  as  those  presented  in  Reference  7. 
This  indicates  the  grids  generated  on  an  IRIS  4D/70GT  computer  system  or 
similar  systems  are  consistent  with  previously  documented  work. 

The  generation  of  the  algebraic  grid  took  863.81  CPU  sec  (30  min 
wall  clock  time).  This  breaks  down  to  1.2  x  10-2  CPU  sec/point.  An  el¬ 
liptic  grid  with  five  iterations  was  also  generated  to  obtain  an  idea  of 
the  time  it  takes  to  generate  an  elliptic  grid.  The  generation  of  the 
elliptic  grid  took  more  than  twice  as  long  as  the  generation  of  the  al¬ 
gebraic  grid.  The  times  are  2225.01  CPU  secs  or (68  minutes  wall  clock), 
which  comes  to  3.75  x  10~^  CPU  sec/iteration/point.  No  elliptic  grids 
are  shown  since  they  demonstrate  little  or  no  difference  from  the  alge¬ 
braic  grids  shown  in  Figures  2a  through  4b. 
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5.1.2  Ogive  -  Cylinder  -  Ogive  with  Fins 

This  configuration  consists  of  an  ogival  nose,  a  cylindrical  raid- 
section  and  an  ogival  boattail  with  fins.  Once  again,  a  sting  extends 
from  the  aft  end  of  the  body.  Figure  5  displays  this  geometry. 

This  configuration  also  uses  a  C-type  system  for  the  field  grid 
around  the  body.  It  consist  of  four  (4)  blocks.  Each  block  has  a  total 
of  33600  (140x24x10)  points,  for  a  sum  total  of  134400  points  for  the 
entire  field  grid.  Figures  6a  through  7b  display  the  side  and  front 
views, of  the  algebraic  grid.  Figure  8  shows  a  perspective  view  of  the 
grid  of  one  of  the  blocks  near  the  body.  Once  again,  the  grids  gener¬ 
ated  on  the  IRIS  4D/70GT  are  comparable  to  those  generated  on  Cray  com¬ 
puter  systems  as  documented  in  References  7  and  18. 

The  generation  of  the  algebraic  and  elliptic  grids  (5  iterations) 
took  2022.25  CPU  secs  (1.5X10-2  CPU  sec/pt)  and  4510.98  CPU  sec 
(3.70x10  3  CPU  sec/iteration/pt) ,  respectively.  The  wall  clock  time  for 
the  generation  of  the  algebraic  grid  was  75  minutes.  No  wall  clock  time 
was  observed  for  the  generation  of  the  elliptic  grid,  although,  it  can 
be  estimated  at  approximately  180  minutes.  Notice  that  the  times  speci¬ 
fied  in  CPU  seconds  per  point  for  the  algebraic  grid  and  CPU  seconds  per 
iteration  per  point  for  the  elliptic  grid  for  the  first  configuration  as 
compared  to  this  configuration  are  similar  although  not  the  same.  The 
difference  can  be  attributed  to  the  larger  number  of  blocks  increasing 
the  1/0  time. 
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5.1.3  Elliptic  Airframe  with  Control  Surfaces 

This  last  configuration  is  of  an  advanced  airframe  consisting  of  an 
elliptical  cross-section  for  the  upper  surface  of  the  body,  a  flat  bot¬ 
tom  surface,  and  vertical  and  horizontal  control  surfaces.  Only  half  of 
this  airframe  was  gridded,  as  seen  in  Figure  9. 

The  field  grid  surrounding  this  airframe  is  a  H-type  system  con¬ 
sisting  of  five  (5)  blocks,  for  a  total  of  111040  points  for  only  half 
of  the  airframe  and  the  surrounding  domain.  The  dimensions  for  each 
block  are  as  follows : 

Block  1  -  70x24x10  (16800)  points 
Block  2  -  70x24x22  (36900)  points 
Block  3  ~  70x16x31  (34720)  points 
Block  4  -  20x24x31  (14880)  points 
Block  5  -  20x16x24  (7680)  points 

Block  1  encompasses  part  of  the  upper  surface  of  the  body  up  to  and 
including  one  side  of  the  vertical  tail.  Block  2  has  the  rest  of  the 
upper  surface  of  the  body  from  the  other  side  of  the  vertical  tail  to 
the  upper  surfaces  of  the  horizontal  control  surfaces.  Block  3  contains 
all  of  the  bottom  surface,  including  the  flat  sides  of  the  horizontal 
control  surfaces.  The  last  two  blocks  -  blocks  4  and  5  -  encompass  the 
entire  domain  in  front  of  the  airframe.  Figures  10a  through  13b  show 
front,  side  and  top  views  of  this  configuration,  and  Figures  13a  and  13b 
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show  a  perspective  view  of  the  entire  domain  and  a  closeup  view  of  the 
grid  near  the  body.  These  grids  are  also  comparable  to  the  grids  docu¬ 
mented  in  References  8  and  12. 

The  algebraic  grid  took  3038.7  CPU  secs  to  generate,  while  the  el¬ 
liptic  grid  (5  iterations)  took  6500  CPU  secs.  With  respect  to  CPU 
time,  the  algebraic  grid  took  2.74x10  ^  CPU  sec/pt  and  the  elliptic  grid 
took  3.8x10  ^  CPU  secs/iteration/pt. 

5.2  Porting 

After  verifying  that  Program  EAGLE  worked  correctly  on  the  IRIS 
4D/70GT  computer  graphics  system,  the  code  was  placed  on  a  SUN  4/280 
computer  system  to  determine  if  the  code  is  portable  across  different 
computer  systems. 

After  placement  of  the  code  on  the  SUN  4/280  computer  system,  the 
FORTRAN77  version  of  the  EAGLE  code  compiled  without  errors  and  without 
having  to  make  any  modifications  or  changes  to  the  code.  This  version 
of  Program  EAGLE  executed  correctly  for  all  three  configurations 
described  earlier.  The  CPU  times  obtained  to  generate  these  three 
configurations  on  the  SUN  4/280  computer  system  are  comparable  to  those 
times  obtained  on  the  IRIS  4D/70GT. 
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Chapter  6 
Conclusions 


The  results  of  this  effort  show  that  the  FORTRAN77  version  of  the 
EAGLE  code  with  its  own  NAMELIST  input  emulator  works  as  well  as  the 
version  written  for  the  Cray  systems.  The  three  configurations  used  to 
test  the  FORTRAN77  version  of  the  EAGLE  code  have  demonstrated  the  capa¬ 
bilities  of  the  code,  as  well  as  obtaining  results  consistent  with  pre¬ 
vious  work  documented  in  References  7  through  9,  12,  and  18.  This 
version  of  the  EAGLE  code  can  then  be  used  on  other  systems  having  the 
FORTRAN77  compiler  with  little  or  no  modifications  to  the  code  as  demon¬ 
strated  on  the  SUN  4/280  computer  system.  As  power  and  speed  of  com¬ 
puter  systems  increase,  the  more  plausible  it  becomes  to  execute  Program 
EAGLE  in  an  interactive  environment,  especially  when  tied  to  a  graphics 
program  allowing  the  user  to  “see"  the  work  being  achieved. 

As  with  any  other  effort,  this  effort  is  not  the  one  to  end  all 
efforts.  But  it  is  a  beginning  of  a  powerful  tool  for  the  computational 
fluid  dynamics  community.  Further  efforts  are  needed  to  give  the  user 
greater  control  and  flexibility  of  the  execution  of  the  code  and  allow 
visualization  by  the  user  of  the  problem  at  hand. 
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f  elliptic  missile  field  grid 


Figure  3b  -  Close-up  front  view  of  elliptic  missile  field  grid 
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Figure  4a  -  Perspective  view  of  overall  elliptic  missile  field  grid 
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Figure  4b  -  Close-up  perspective  view  of  elliptic  missile  field  grid 


Figure  5  -  Ogive-Cylinder-Ogive  with  fins  surface 
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Figure  6b  -  Close-up  side  view  of  Ogive-Cylinder-Ogive  with 


Figure  7a  -  Front  view  of  Ogive-Cylinder-Ogive  with  fins  field  grid 
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Figure  7b  -  Close-up  front  view  of  Ogive-Cylinder-Ogive  with  fins  field  grid 


Figure  9  -  Elliptic  airframe  with  control  surfaces  surface  grid 
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Figure  10a  -  Front  view  of  elliptic  airframe  field  grid 
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ii 


Figure  10b  -  Close-up  front  view  of  elliptic  airframe  field  grid 
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Figure  13b  -  Close-up  perspective  view  of  elliptic  airframe  overall  field  grid 
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Appendix  A 

SAMPLE  PROGRAM  SET-UP  LISTING 
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nno  ooooon  non  non 


L> 


c 

c 

c 

c 

c 


CHARACTERS  NAMLST  >  NAML.XS  ;  CD  A  TA 
CHARACTER# t 44  STRING 
CHARACTEE#64  ICHR 

INTEGER  CPGS  t  IPGS  >  RPGS  >  PRMPT 
INTEGER  UP  ?  DGWN  >  RSTAR7 

INTEGER  PI  INFS  >  PINL  ?  PITEXT  >  PIECHO  t  PLS  >  Pf.T  ;  PHI.  ?  PNJ  > 
i  PIcQUAL  >  PNUM 

PARAMETER (  Nf.HAR  =  31  y  NINTG  =  27  f  NREAL.S  -  Vt  > 

PARAMETER <  NCVAL  =  36  >  NIVAL  =  2#NVALMX+  4 ,1  >  NKVAI.  -  »t«P3  t  5  ) 
PARAMETER  (  NX  HAT  =  3  >  NEDA7  =  ,?  ) 

PARAMETER (  NVAR  =  NCHAR  +  NINTH  +  NREALS  ) 

PARAMETER (  NJ5AT  =  NIDAT  t  NRDAT  ) 

PARAMETER (  IIP1  =  NVALMX  +  41  >  ICP1  =  2  ) 


PARAMETER  (  RNONE  =  3*402824E+3o  ,  INDNE  =  2147483647  ) 


DIMENSION  NAMLIS (  NVAR  )  >  NPIMJ  (  NVAR  )  ,  NDIM,?(  NVAR  >  t 
S  NAMLST (  2  ) 

DIMENSION  ICHR <  NCVAL  )  ?  INTGR<  NIVAL  )  y  REAL  SC  NRVAL  ) 
DIMENSION  CDATAC  NDAT  )  f  JDATC  NIDAT  )  ?  RDATC  NRDAT  ) 
DIMENSION  CPOSC  NCHAR  )  t  IPGSC  NINTH  )  >  RPOS'.  NREALS  ) 


THIS  SET  FOR  THF  CHARACTER  V ART ARLES 


EQUIVALENCE 

< 

ICHR  <  1) 

ITEM 

) 

f 

( 

ICHR (  2>  > 

CLASS 

) 

EQUIVALENCE 

( 

ICHR (  3) 

*  INTERP 

) 

f 

( 

ICHR (  4)  y 

PROTYP 

) 

EQUIVALENCE 

( 

ICHR <  5) 

FUN 

) 

y 

( 

ICHR (  6)  > 

BLEND 

) 

EQUIVALENCE 

( 

ICHR (  y> 

REWIND 

> 

y 

( 

ICHR ( 10 )  > 

CGN7YP 

) 

EQUIVALENCE 

< 

XCHRCn  ) 

RESTART 

) 

y 

( 

ICHR<J2)  y 

CHECK 

) 

EQUIVALENCE 

< 

ICHR (13) 

SMOCON 

) 

y 

< 

I CHR (14)  y 

ALL 

) 

EQUIVALENCE 

C 

ICHR<15) 

WEIGHT 

) 

y 

( 

ICHR ( 16 )  y 

MATH 

) 

EQUIVALENCE 

C 

ICHRC17) 

OUTER 

) 

y 

( 

ICHR (18)  / 

OFIRST 

) 

EQUIVALENCE 

< 

ICHR ( 3  9 ) 

PCROSS 

) 

y 

( 

J  CHR ( 20 )  y 

PROFOE  t  PROF 

EQUIVALENCE 

( 

ICHR (23) 

CONCIJT 

> 

r 

( 

ICHR (24)  f 

CONTIN 

) 

EQUIVALENCE 

( 

ICHR (23) 

TRIAD 

) 

y 

( 

ICHR  <  26 )  f 

NEW 

) 

EQUIVALENCE 

( 

ICHR (27) 

KSTORE 

) 

y 

( 

ICHR (28)  > 

FORM 

) 

EQUIVALENCE 

c 

ICHR ( 29 ) 

MODE 

) 

y 

( 

ICHR (30)  y 

I PRINT 

) 

EQUIVALENCE 

( 

I CHR (31) 

PART 

) 

y 

( 

I CHR ( 32 )  > 

RLKERR 

) 

EQUIVALENCE 

( 

ICHR (33) 

INTORT 

) 

t 

( 

ICHR (34)  , 

F II.  NAM 

) 

EQUIVALENCE 

( 

ICHR (33)  * 

CONTUR 

) 

THE  NFXT  SET  IS  FOR  THE  INTEGER  VARIABLES* 


EQUIVALENCE  <  JNTGRC  1)  #  START  >  POINT  ) 
EQUIVALENCE  C  I NT HR (  4)  y  ISTART  >  IPOTNT  > 


64 


EQUIVALENCE 

( 

INTGRC  7) 

7 

END 

*  <  IN7GR< 10  J 

7 

I  END 

5 

EQUIVALENCE 

( 

INTGRUS) 

t 

BLOCK 

B  7  BLOC  ) 

EOUI VALENCE 

< 

INTGR<  3  4 ) 

7 

r  BLOCK 

IB  ) 

EQUIVALENCE 

( 

IN  TOR (IS) 

7 

ORDER 

SMODIR  7  S  ) 

EOUSVAI  FNC" 

< 

JNTGRUS) 

t 

RORDER 

SMOCOM  7  RS  ) 

EQUIVALENCE 

C 

INTGR  <?  t ) 

9 

OPOINT 

) 

9  <  INTGR  <  74 ) 

7 

UMAX 

) 

EQUIVALENCE 

V 

INTGR  <7.5 ) 

7 

FILE 

) 

7  <  INTGR ( ?o ) 

7 

CONUP.T 

> 

EQUIVALENCE 

( 

INTGR (77  > 

VALOUT • VALUE  )  .  <  INTGR < 

73 )  7  LOCAT 

EQUIVALENCE 

< 

INTGR ( 31 ) 

* 

SEGMENT 

) 

*  (  INTGR  <  32  > 

7 

POINTS 

> 

EQUIVALENCE 

< 

INTGR ( 33 ) 

7 

DIRECT 

) 

f  <  INTGR  <  .14 ) 

7 

NDEX 

) 

EQUIVALENCE 

< 

INTGR  <  35) 

7 

PROn.IR 

9 

PRuJ  ) 

EQUIVALENCE 

< 

tNToRCSS) 

9 

SIZE 

) 

EQUIVALENCE 

( 

JNTGRC4I ) 

7 

TERMS 

) 

7  <  INTGR < I  IF I 

t  X  7  E  R  il  S 

) 

THE  Nr 

XT  SFT  IS  F 

OR 

THE  REAL 

VAR J ABLER* 

EQUIVALENCE 

< 

REALS <  I) 

9 

ACCPAR 

) 

7  (  REAL  S  <  7 ) 

7 

TOL 

) 

EQUIVALENCE 

( 

REALS <  3> 

9 

CONEAC 

) 

7  <  REALS <  4) 

f 

SPAVAL 

) 

EQUIVALENCE 

( 

REALS <  5) 

9 

VALUES 

) 

c 

C  THIS  SET  OF  EQUIVALENCE  STATEMENTS  MAINLY  USED  WHEN  CHARACTER 

C  VALUES  ARE  INPUTS  TO  INTEGER  OR  REAL  VARIABLES  * 

C 

EQUI VALENCE (  IDATU)  t  UP  )  7  <  IDATC2)  t  DOWN  ) 

EQUIVALENCE (  IDATCS)  r  RSTART  ) 

C 

EQUIVALENCE <  RDATU?  7  OPTIMUM  )  t  <  RDATC2)  7  EXTRAP  ) 


C 

c, 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


THIS  SFT  OF  STATEMENTS  ARE  FOR  QUANTITIES  OF  TYPF  CHARACTER 


DATA  NAMLST/  'GINPUT'  7  'GOUTPUT'  / 


DATA 

(NAMLISU)  f  J»J 

jNCHAR)/  'ITEM 

'7  'CLASS'  7 

'INTER?' 

7 

£ 

' PROTYP ' 

7 

'FUN'  f 

'BLEND'  f 

'REWIND' 

7 

% 

'CONTYP ' 

7 

'RESTART'  7 

'CHECK'  7 

' RMOCON ' 

7 

s 

9  ALL 9 

t 

'WEIGHT'  7 

'  MA7H '  7 

'OUTER' 

7 

% 

9 DEIRST ' 

f 

' DCROSS '  j 

'PROFOR'  7 

'CONCUT' 

7 

i 

' CGNTIN 9 

7 

'TRIAD'  7 

'NEW '  7 

' KSTuRE ' 

7 

s 

9  FORM ' 

7 

'MODE'  f 

'IPRINT'  7 

'PART' 

7 

s 

9 3LKERR 7 

9 

'INTORT'  7 

'FILNAM'  7 

'CONTUR' 

/ 

THIS  SET  OF  STATEMENTS  ARE  FOR  QUANTITIES  OF  TYPE  INTEGER 


DATA  (NAMI.  IS<  I )  t  I-NCHAR+X  7 NINTG+NCHAR )  / 


& 

'START' 

7 

'ISTART' 

7 

'END'  7 

' IEND ' 

7 

2 

'BLOCK' 

7 

'  IBLOCK ' 

7 

'ORDER'  7 

'SMODIR' 

7 

% 

'RORDER' 

7 

'SMOCOM' 

7 

'POINT'  7 

' IPOINT ' 

7 

s 

'OPOINT' 

7 

'  ITMAX ' 

7 

'FILE'  7 

' CGNUPI ' 

7 

s 

'VALOUT' 

7 

'LOCAf ' 

7 

'SEGMENT'  7 

'POINTS' 

7 

& 

'DIRECT' 

7 

'NDEX' 

7 

'  PROD IR '  7 

'SIZE' 

7 

s 

'TERMS' 

7 

' I TERMS' 

7 

'VALUE'  / 

THI 

3  SFT  OF  TSTA7 

EMENTS  ARE  FOR 

QUANTITIES  GF 

TYPE  RE 

AL 

DATA 

(NAMLISU)  7l=N 

iNTG*NCHAR+J 

»hi:har+nintgtNreas„s>  / 

a 

'ACCPAR' 

V 

'TOL' 

7 

'  CONFAC '  7 

'SPAVAL' 

7 

s 

'VALUES' 

/ 

C 

C  THIS  FIRST  SFT  OF  DATA  VALUES  FOR  THE  ARRAY  NRTM  IS  FOR  CHARACTER 
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ooo  n  n  n  onnoo  o  o  ooo  o  o  o  o  o  n  o  c;  o  n  o  o  o  o  o  o 


VARIABLES .  ANY  ADDITIONS  Or  CHARACTER  VARIABLES  SHOULD  COME  AT  THE 
END.  IN  THF  EVENT  AN  ARRAY  OF  VARIABLE  SITE  IS  TO  BE  ADDED  MAKE 
SURE  THOSE  ARRAY  SIZES  ARE  ALWAYS  AT  THE  END. 

THE  SECOND  SET  IS  FOR  INTEGERS.  NOTE  THAT  THE  VARIABLE  DIMENSIONED 
ARRAY  SIZES  ARE  AT  THE  END  OF  THE  INTEGER  LIST. 


THE  LAST  GRCUF  Or  VALUF3  ARE  FOR  REAL  VARIABLES.  AGAINt  NOTE  THAT 
THE  VARIABLE  DIMENSIONED  ARRAY  SIZES  ARE  AT  THE  END  OF  THE  REAL 
LIST. 


DATA  NDIM.I/  It  If  It  It  It  it  Xt  it  it  it  it  it  i  t  it  it  it  it  it 
S  b  If  tr  h  If  h  If  b  li  h  It  li  li 


It  It  It  It  It  It  It  it  If  If  it  It  If  If  if  if 
it  it  if  it  it  if  it  if  it  it  it 


£ 


If  It  If  if  1  / 


DATA  NDIM7/  It  Jt  It  It  It  3t  If  it  it  it  It  It  it  It  It  If  if  3t 
£  it  i»  it  it  it  if  it  it  if  it  it  if  I or 1 t 


£  3t  3t  3t  3t  It  If  3t  3f  3t  3t  3t  3f  If  If  if  it 

S  It  It  li  It  tf  It  It  it  NVALIiX t NVAI.MX t  it 


£ 


It  It  It  i t  DIMP3/ 


THESE  ARRAYS  USED  BY  THE  PARSER  FOR  POINTING  TO  A  MEMORY  LOCATION. 


DATA  CPOS/  it  ?.t  3 1  4 1  5t  6t  9 1 10/11 1 12/  J  3/  1 4 1 15f  1 6/  17t  i  8 1  i9t  20 1 
£  23/24/25/26/27/28/29/30/31/32/33/34/35/ 

DATA  IPOS/  It  4 1  7t  10tl3f  14fl5fi5flSfl8f  it  4/7.1/74/75/26/27/28/ 
£  31 / 32/ 33/34 /35 / 38/ 4 1 / I IP1 t 27/ 

DATA  RPOS/  it  2f  3t  4t  5/ 


THESE  STATEMENTS  ARE  ADDITIONS.  070988  <*U*  AM> 

DATA  IIP/  2147483644  /  t  DOWN/  2147483645  / 

DATA  RSTART/  2147483646  / 

DATA  0PTIMIIM/2147483647.1/  f  EXTRAP/2147483647 ..1/ 

DATA  CDATA/  'UP'  t  'DOWN'  t  'RESTART'  t  'OPTIMUM'  t  'EXTRAP'  / 


DATA  PLINES/  1  /  t  PIECHO/  0  /  t  PINL/  1  /  t  Pl.S/  1  /  t  PL  7/  1  /  t 
£  PNI./  1  /  t  PN.I/  1  / 

C 

C 
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n  n  n  n  o  o  o  o  o  o  o  n  o 


n 

RFAA  NAHFLlST  *  5  ************  **  f  ?  5  *  M  S*  *  3  25  5  I*  *  ** 
C 

C  THE  FOLLOWING  LINES  REPLACE  THE  LINE  *  READ  INPUT-*  THESE  LINES  CALL 
C  THE  NAMELIST  INPUT  EMULATOR .  MOTE  THAT  THE  ROUTINE  "  PROMPT  * 

C  SETS  UP  A  PROMPT  FOR  THE  INTERACTIVE  USER  * 

r 

PRMPT  =  PRMPT  +  i 
MAS  =  0 
C 

IF (  IBATCH  .EG.  0  >  CALL  PROMPT*  PRilPT  ) 

NOW  SFAfl  IN  THE  INPUT  STRING  ANA  An  TEEii  I  WE  IF  ITS  TOO  LONG. 

1—  CONTINUE 

CALL  R ASTRO*  STRING  7  LOT  »  RFAAFL.  r  NSTRNG  ?  XSTRNO  :  IS  TAT  ) 

IF'  ISTAT  .EG.  X  )  GO  TO  300 

PASS  THE  CHARATER  STRING  ‘STRING  *  70  7Hn  PARSER  ANA  PARSE  THE  INPUTS. 

CALL  PARSER*  STRING  r HAUL  ST  <  3 ) vNAHLISrNDXKXf NSIKSfr XCHRr XNIGRr 

1  .REALS ; GRATA* I  RAT  7 RAAT  ?  CPOS  7  IPOS  5  RPOS  7  NVAR  7 

%  NCHAR  t  NINTG  t  NREALS  t  NCVAI-  t  NIVAL  ?  NRVA1  7 

S  NAA  f  ?  MI  DAT  ?  NRAAT  7  NAS  7  PLXNE8  t  PINl.  7 

2  PITEXT  t  PIECHO  7  PLS  t  PLT  ;  PNL  ?  PNJ  7  PIEGUAL  > 

2  PNIJM  7  ISTAT  ) 

IF  *  ISTAT  .EG.  3  )  GO  TO  100 


IF <  NAS  .EG.  0  )  GO  TO  1?? 


67 


Appendix  B 

GENERALIZED  PARSER/NAMELIST  INPUT  EMULATOR  LISTING 
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SUBROUTINE  RDSTRG (  STRING  >  LCT  t  IREAD  t  NSTRNG  ?  JSTRNG  r  I ST AT ) 
CHARACTERS < * )  STRING 
C 

IF<  IREAD  ♦ EQ  *  1  )  THEN 
READ (  2r?9 V  )  STRING 
ELSE 

READ (  5>?9?  )  STRING 
ENDIF 
C 

LCT  =  ILEN <  STRING  >  ISTRNG  ) 

C 

IF (  LCT  *  GT  t  NSTRNG  )  THEN 
LST  =  LCT  -  NSTRNG 

PRINT  *t'  ERROR  -  TOO  MANY  CHARACTERS  IN  THIS  LINE  ' 

PRINT  *t  STRING (  l t  LCT  ) 

PRINT  %i'  SHORTEN  LINE  BY  APPROXIMATE  Y  ' *LST» '  CHARACTERS ♦ 

PRINT  *9'  ' 

ISTAT  =  1 
RETURN 
ENOIF 
C 

999  FORMAT  <  A  > 

C 

RETURN 

END 
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SUBROUTINE  PARSERC  STRINGrNAMLSTt NAHLIS»NflUiJ  f NDXMJ?fXCHR  f  XMTGRr 
l  REALS  *  Oft AT  A  t  TDAT  f  RDAT  t  CPOS  t  IPOS  t  RPOS? 

J  NVAR  >  NCHAR  »  NINTH  >  NRFALS  ?  NOVAI  >  NJVAI.  f 

z  NRUAL  r  NDAT  >  Nt  DAT  >  NRDAT  >  NDS  >  L  >  INI.  > 

l  ITEXT  t  IECHO  f  I.S  ?  I. T  t  Nl.  ?  NJ  t  XEQUAI.  > 

3  NUM  f  TSTAT  ) 

c t#***********^***^***^************************************************1 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


'HESE  STATEMENTS  USED  FOR  PARSING  ♦  0709BB  <*(  is  AM>  CiADCTljri 

THE  FOLLOWING  IS  a  DESCRIPTION  OF  THE  UARIABI.es  USr.D  FOR  PARSING. 

STRING  -  INPUT  CAEACTFR  I.  INF t  PASSED  TO  THE  PARSER  FOR 
PARSING  * 

NAHLST  -  CONTAINS  THE  NAME  OF  THF  NAMELIST. 

NAML IS  -  CHARACTER  ARRAY  DIMENSIONED  TO  THE  NUMBER  OF 
VARIABLES  IN  THE  NAMELIST  (  NVAR  )  AND  WHICH 
CONTAINS  THE  LIST  OF  INPUT  ITEMS  IN  NAMELIST. 

THIS  ARRAY  IS  ORGANIZED  WITH  INPUT  ITEMS  OF  TYPE 
CHARACTER  FIRST  >  INTEGER  SECOND  >  AND  REALS  LAST. 

ICHR  -  CHARACTER  ARRAY  WHICH  HOLDS  THE  VALUES  FOR  THE 

INPUT  ITEMS  OF  TYPE  CHARACTER*  THIS  ARRAY  IS  THEN 
EOUIVALENCED  TO  THE  APPROPRIATE  VARIABLES.  THIS 
ARRAY  IS  DIMENSIONED  TO  THE  NUMBER  OF  CHARACTER 
VALUES  <  NCVAt.  )  ♦ 

INTGR  -  INTEGER  ARRAY  WHICH  HOLDS  INTEGER  VALUES  FOR  THE 
INPUT  ITEMS  OF  TYPE  INTEGER .  THIS  ARRAY  IS  THEN 
EQUIVALENCE  TO  THE  PROPER  INTEGER  VARIABLES* 

INTGR  IS  DIMENSIONED  TO  THE  NUMBER  OF  INTEGER 
VALUES  (  NIVAL.  )  IN  NAMLIS. 

REALS  -  REAL  ARRAY  WHICH  HOLDS  REAL  VALUES  FOR  THE  INHUT 

ITEMS  OF  TYPE  REAL .  IT  THEN  IS  EQUIVALENCE!!  TO  THE 
APPROPRIATE  VARIABLES.  THE  ARRAY  REAL  IS  DIMENSIONED 
TO  THE  NUMBER  OF  REAL  VALUES  <  NRVAL  )  IN  NAMLIS. 

NDIM1  t  NDIM2  -  ARRAYS  CONTAINING  THE  LOWER  X  UPPER  BOUNDS  » 

RESPECTIVELY >  OF  EACH  OF  THE  VARIABLES  IN  NAMLIS  IN 
THE  ORDER  OF  CHARACTER  FIRST*  INTEGER  SECOND  >  AND 
REAL  AS  LAST.  NDIMt  AND  N0IM2  ARE  DIMENSIONED  BY  THE 
THE  TOTAL  NUMBER  OF  VARIABLES  IN  NAMELIS<  NVAR  ). 

NCHAR  -  NUMBER  OF  CHARACTER  VARIABLES 

NINTG  -  NUMBER  OF  INTEGER  VARIABLES 
NRFALS  -  NUMBER  GF  RFAL  VARIABLES 

»;>NOT£:  NVAR  =  NCHAR  +  NINTG  +  NREALS<<<<<<< 

THE  FOLLOWING  ARRAYS  ARE  USED  WHEN  CHARACTER  ^**FORMATXOH  X8  PASSED 
IN  INTEGER  or  REAL  VARTARLES<  IN  THE  INPUT  >  INSTEAD  ^INTEGER 
nfc  t- l- a l  VAUJFS.  THIS  IS  USEFUL  IN  THAT  IT  ALLOWS  THE  USFR  SOME 

FI  EXXBU  TTYf  A  .e.f  THE  USER  DOES  NOT  HAVE  TO  <5UPPi:JvTt,!TpeT 
PR I ATE  VALUES  (  INTEGER  OR  REAL  )  BUT  CAN  INPUT  SAY  FIRST  > 
^AST'/oR  'SAME'  <  OR  WHATEVER*  AS  DESCRIBED  BY  THE  USERS  MANUAL 
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C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


’PROGRAM  EAGLE  -  USERS  MANUAL  VOLUME  II*  AY  DR*  JOE  F*  THOMPSON 
/IT SO*  STATE  UNIV  t  OEPT  OF  AEROSPACE  ENGINEERING  r  MS  >  * 


CDATA  -  CHARACTER  ARRAY  CONTAINING  THE  NAMES  OF  CONSTANTS. 
THIS  ARRAY  IS  DIMENSIONED  TO  THE  NUMBER  OF 
INTEGER  AND  REAL  CONSTANTS  <  NDAT  >. 

IDAT  -  INTEGER  ARRAY  CONTAINING  INTEGER  CONSTANTS 

INITIALIZED  BY  THE  MAIN  PROGRAM  *  THE  ARRAY  IS 
EflUIVAt.ENCED  TO  THE  APPROPRIATE  VALUES*  IDhT  IS 
DIMENSIONED  BY  THE  NUMBER  OF  INTEGER  VALUES  IN 
*  COMMON/  CONST  /  *.**<  NJ  DAT  )* 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 


RDA7  - 


»»»>NOTS 


REAL  ARRAY  CONTAINING  INTEGER  CONSTANTS  INITIALIZED 
IN  THE  MAIN  PROGRAM*  THIS  ARRAY  IS  EQUIVALENCES  TG 
THE  APPROPRAITE  VALUES  AND  IS  DIMENSIONED  TO  THE 
NUMBER  OF  REAL  VALUES  IN  'COMMON/  CONST  / 

(  NRDAT  )  ♦ 


CPOS  -  INTEGER  ARRAY  WHICH  HOLDS  THE  ACTUAL  ARRAY  POSITIONS 
OF  ALL  THE  NAMELIST  INPUT  VARIABLES  OF  TYPE 
CHARACTER.  CPOS  ACTS  AS  A  POINTER  FOR  THE  PARSER* 


C 

C 

C 

c 

c 

c 

c 

c 


IPOS  -  INTEGER  ARRAY  WHICH  HOLDS  THE  At  ,AL  ARRAY  POSITIONS 
OF  ALL  THE  NAMELIST  INPUT  VARIABLES  OF  TYPE  INTEGER* 
IPOS  ACTS  AS  A  POINTER  FOR  THE  PARSER. 

RPOS  -  INTEGER  ARRAY  WHICH  HOLDS  THE  ACTUAL  ARRAY  POSITIONS 
OF  ALL  THE  NAMELIST  INPUT  VARIABLES  OF  TYPE  REAL  * 
RPOS  ACTS  AS  A  POINTER  FOR  THE  PARSER* 


OTHER  additions: 


C 
C 
C 

c 
c 
c 
c 


ISTAT  -  STATUS  FLAG.  MSTAT  =  O’  NO  ERRORS?  * ISTAT  =  X  * 
INDICATES  AN  ERROR.  IT  WILL  PROMPT  THE  USER  IN 
INTERACTIVE  MODE . OTHERWISE f  EXECUTION  WILL  STOP. 


#*#*##########**#####*######*##*##*##*#***#####**###***#######**####### 


CHARACTER#  <40  NAML.ST  ?  NAMLIS  i  CDATA 
CHARACTER#t44  TEXT 
CHARACTER# < * )  STRING 
CHARACTER# ( * )  ICHR 
C 

INTEGER  CPOS  r  IPOS  t  RPOS  t  POINTR 
C 

PARAMETER (  NT EXT  =  J 44  ) 

C 

DIMENSION  I CHR <  NCVAL  )  r  INTGR (  NJVAl  >  ?  REALS (  NKVAL  ) 
DIMENSION  NAMLIS (  NVAR  )  /  NDIMl (  NVAR  )  f  N0IM2C  NVAR  ) 
DIMENSION  CDATA (  NDAT  )  r  IDAT <  NIDAT  )  f  RDAT (  NRDAT  ) 
Dli.'NSION  CPOS<  NCHAR  )  ,  IPOSC  NINTH  )  >  RPOSC  NREALS  ) 

C 

C*##**#**>r*  v#  ###**### 
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c 

ICO  =  ICHAR f  'O'  ) 

IC?  =  ICHAR  <  ) 

ICA  -  ICHAR (  #A'  ) 

.CC7:  =  I  CHAR  (  > 

ICQU  =  ICHAR  (  ""  > 
ICCM  =  ICHAR (  ) 

ICFOS  =  1CHAR<  ) 

1  CM IN  =  ICHARC  '-'  ) 
ICPER  =  ICHAR  <  ) 

C 

IP AREN  «  0 
IROINT  =  1 
C 

LNL  -  Il.EN<  NAML3T  a  ) 
m.  =  LNL  +  2 


DO  1  I  =  ImNTEXT 

i  text <  x:x  )  =  '  ' 

C 

ISTRNG  =  LEN <  STRING  ) 

LC  =  ILEN<  STRING  t  TSTRNG  ) 

C 

IF  <  STRING <  1U  )  ♦  NE ♦  '  '  )  THEN 
C 

IF (  STRING (  Ui  )  . EG  ♦  'C'  )  THEN 
C 

PRINT  *?  STRING <  1JI.C  ) 

PRINT  tf '  ' 

GO  TO  999? 

C 

ELSE  i.F<  STRINGS  111  )  » EQ#  ' E '  •  AND»  I.  *En*  X  )  THPN 
C 

PRINT  Xt  STRING <  1  JLC  ) 

PRINT  %t  9  9 
IECHO  =  1 
C 

ELSE  IF <  L  ♦ EQ ♦  1  )  THEN 
C 

PRINT  %t  '  ERROR  -  COLUMN  1  MUST  EITHER  BE  BLANK*  HAVE  A  ' f 
%  '<C>  OR  AN  <E>  ' 

PRINT  *t  '  9 
ISTAT  =  1 
GO  TO  9999 
C 

ELSE 

C 

PRINT  Xt9  ERROR  -  COLUMN  1  MUST  BF  BLANK  ' 

PRINT  Xf '  ' 

ISTAT  =  1 
GO  TO  9999 
C 

END  IF 
C 

ELSE  IF <  STRING (  1 ; 1  )  ♦ EQ ♦  '  '  ♦AND,  IECHO  ♦EO,  X  )  THEN 
C 

PRINT  *>  STRING (  1JLC  > 

PRINT  Xf '  ' 

C 

ENDIF 
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IF (  L  *  EG*  1  )  THEN 

IF <  STRINGC  2J2  )  *NF ♦  )  THEN 

PRINT  %?'  ERROR  -  DHL. LAP  SIGN  REQUIRED  IN  COLUMN  2  ' 

PRINT  *>  '  ' 

IS  TAT  =  1 
60  TO  9999 
ENDIF 

NDS  =  INDEX  (  STRING  C  3:LC  )>'$') 

IF  <  NDS  « NF  *  0  )  NDS  =  NDS  +  2 
ELSE 

NDS  =  INDEX  (  STRING  (  2U.C  )  r  9%9  > 

IF (  NDS  *NE*  0  )  NDS  =  NDS  +  1 

ENDIF 

IF <  NDS  ♦ EG ♦  0  )  THEN 
LS  »  LC 
ELSE 

LS  =  NDS 
ENDIF 

TEST  TO  SEE  IF  THE  PROPER  NAME  Or  THE  NAMELIST  IS  FOUND  IN  THE 
STRING  IN  THE  RIGHT  PLACE  * 

IF  C  L  ♦ EG  *  I  )  THEN 

IF<  STRING (  3 1 NL  )  * NE ♦  NAMLST  )  THEN 

PRINT  ERROP  -  ' t NAMLST t '  NOT  FOUND  BETWEEN  COLUMNS  o  THRU' 

S  NL  * '  ♦  '  j  STRING <  KNL  ) 

PRINT  ' 

ISTAT  =  1 
GO  TO  9999 
ENDIF 

TEST  TO  SEE  IF  A  BLANK  EXISTS  AFTER  NAMLST  IN  STRING ♦ 

IF<  STRING <  NL+ltNL+J  )  *NE*  '  '  )  THEN 

PRINT  %$'  ERROR  -  COLUMN  AFTER  9 1 NAMLST r  '  MUST  BE  BLANK  ♦'> 
l  STRING <  UNt.+  X  ) 

PRINT  %f  '  ' 

ISTAT  =  1 
GO  TO  9999 
ENDIF 

SET  STRING  POINTER (  NL  )  TO  PROPER  POSITION  FOR  THE  PARSER  . 

NL  =  NL  +  2 
ELSE 
NL  -  2 
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C 

END  IF 
C 

10  CONTINUE 


THIS  SECTION  PUTS  TOGETHER  A  STRING  OF  CHARACTERS  TO  FORM  A 
QUANTITY  TO  BE  LATER  DETERMINED  IF  THE  QUANTITY  IS  IN  THE 
LIST  NAML  J3 ♦ 

DO  20  I  =  NL  t LS 


TEST  TO  SEE  JF  ANY  LEADING  BLANKS  EXIST t IF  SO  INCREMENT 

IF (  STR JNG  <  I J I  )  ♦ £3  *  '  '  ♦AND*  TEXT  <  } : 1  )  « EG  *  '  '  >  GO  TO  20 

IF  <  STRJNG<J:i)  ♦ EG ♦  ♦ AND  *  J  ♦ FQ ♦  LS  >  THFN 

IF <  TEXTUJJ)  *NE*  '  7  )  THEN 

PRINT  ft'  ERROR  -  NOT  FINISH  READING  THE  QUANTITY  WHEN  '> 

X  'INPUT  LINE  WAS  TERMINATED  ' » STRING (  L3-7U.S  ) 

ISThT  =  1 
GO  TO  9999 
ELSE 

GO  TO  9999 
ENDIF 

ENDIF 


ICHS  *  ICHAR <  STRING!  III  )  ) 

TEST  TO  SEE  IF  THE  CHARACTER  STRING  STARTS  WITH  AN  INTEGER t  IK  SO 
THEN  ITS  AN  ERROR, 


IF C  L  ♦ EQ ♦  1  >  THEN 

IF  <  ( ICHS  .  GE ,  ICO  .AND,  ICHS  ,  I.  E  »  IC9)  ,  AND  .  TEXT  ( 3  1 1 )  .  EQ  .  '  ') 
X  THEN 

PRINT  # r ' ERROR  -  QUANTITY  CANNOT  START  WITH  A  NUMBER  't 
X  STRING (  1 : 1+7  ) 

PRINT  *r '  ' 

ISTAT  =  1 
GO  TO  9999 
ENDIF 


ELSE  IF <  < (  ICHS  ♦ GE ♦  ICO  .AND,  ICHS  , LE ♦  IC9  )  ,0R, 

X  ICHS  ,  EQ ,  ICGIJ  .OR,  ICHS  » EQ .  ICCM  .OR, 

2  ICHS  « EQ  *  ICPOS  » OR ♦  ICHS  , EQ ♦  ICMIN  .OR. 

X  ICHS  .EG.  ICPER  )  .AND,  TEXT'ltl)  ♦ EQ ♦  '  '  )  THEN 

IEQUAL  =  1 
NL  =  l 
GO  TO  210 

ENDIF 


FIRST  TEST  TO  SEE  JF  THE  CHARACTER  STRING C  IJI  /  IS  BETWEEN  0->9 

AND  A->7.  IF  SO  CONCATENATE.  IF  NOT  THEN  TEST  TO  SEE  IF  AN  EMBBEDDED 
OR  A  TRAILING  BLANK  EXISTS. IF  SO  CONCATONATE.  IF  NOT  THEN  TEST  TO 
SEE  IF  THF  QUANTITY  STARTS  WITH  A  SYMBOL , IF  SO  -  ERROR,  IF  NOT 
THEN  DETERMINE  JF  ANY  OHARACTFR  FOLLOWS  THF  QUANTITY  TO  TERMINATE 
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SEARCH  AND  CONCATENATION 


IF  (  (  ICHS  *  OF  •  ICO  *  AND  *  ICHS  *I.E*  JC?  )  ♦OR* 

(  ICHS  ♦ OF  *  ICA  ♦ AND  *  ICHS  *  LE  ♦  IC7.  )  )  THEN 

lr<  TFXT (  3:i  )  *cQ*  '  '  )  THEN 
TEXT (  1M  )  =  STRING (  IJI  ) 

IL  =  1 
ELSE 

TEXT  =  TEXT (  J  t Jt.  )//STRING(  IJI  ) 

IL  =  It.  +  1 
ENDIF 

ELSE  IF  (  STRING  (  IJI  )  ♦  EH .  '  '  .AND,  TEXT  (  UX  >  .Nr*  ' 


C 


C 

c 


c 

c 


c 

c 


c 

c 

c 


TEXT  ^  TEXT  (  XI 1 1.  ) //STRING  <  IJI  ) 

IL  =  It.  +  1 

ELSE  IF (  STRING (  IJI  )  ♦ EQ ♦  ♦ OR ♦  STRING <  JIT  )  . EG *  ' 

2 

IF (  XPAREN  .EQ.  0  )  THEN 

IF (  STRING (  J : X  )  ♦ EG ♦  '<'  ) THFN 
TEXT  =  TEXT (  1?IL  )//STRING(  I*t  ) 

IL  =  IL  +  1 
IPAREN  -  XPAREN  +  1 

ELSE  IF <  STRING*  IJJ  )  *EG.  ')'  >  THEN 

PRINT  ERROR  -  HISSING  LEFT  PARAENTHESIS  7  > 

%  STRING  < 1-9 : 1 ) 

PRINT  tf  '  ' 

ISTAT  =  1 
GO  TO  9999 
END  IF 

ELSE  IF C  IPAREN  *  EH  ♦  .t  )  THEN 

IF (  STRING (  I X I  )  >  THEN 

PRINT  tf'  ERROR  -  HISSING  RIGHT  PARENTHESIS  'f 
l  STRING (  1-9  J I  ) 

PRINT  tt'  ' 

ISTAT  =  1  > 

GO  TO  999V 

ELSE  IF (  STRING (  IJI  )  * EO ♦  ) THEN 

TEXT  =•*  TEXT <  1JIL  )//STRING<  I J I  ) 

IL  =  IL  +  1 
IPAREN  =  IPAREN  +  1 
ENDIF 

ELSE 

PRINT  tt '  ERROR  -  POSSIBLY  TOO  HANY  PARENTHESIS  '> 
X  STRING \  1-9: I) 

PRINT  tf  '  ' 

ISTAT  =  I 
GO  TO  9999 

ENOTF 

ELSE  IF (  STRING (  I J  J  )  .NE.  '  '  .AND.  TFXT (  ill  )  . EO .  ' 


' )  THEN 


)'  ) 
THEN 


')  THEN 
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STRING <  IJI  )  ♦  NE.  )  THEN 

PRINT  *f'  ERROR  -  QUANTITY  STARTING  WITH  A  SYMBOL  '/ 

3  STRING (  i:i+7  ) 

PRINT  *f'  ' 

ISTAT  =  1 
GO  TO  9999 
ENPIF 
C 

ELSE 

IF  <  STRING  <  I  :i  )  ♦  NE  ♦  '  '»AND.  TEXTUJl)  ♦  NE  ♦  '  ')  GO  TO  SO 
ENPIF 
C 

20  CONTINUE 
C 

IF<  ND3  *  EG  *  0  )  THEM 
L  *  L  +  1 
GO  TO  9999 
Ex*OIF 
C 

30  ;*•  '  NUE 


SET  PC  I?*  "ER  TO  I.  AST  POSITION  WHICH  IS  WHERE  THE  NEXT  CHARACTER 
BEGINS <  MOST  LIKELY  AN  EQUAL  SIGN  )* 


NL  =  I 

LT  =  ILEN<  TEXT » NTEXT  ) 

AFTER  FINDING  PARENTHESIS  THIS  INDICATES  THAT  AN  ARRAY  ELEMENT  IS 
BEING  SPECIFIED  *  THEREFORE  f CONVERT  THE  CHARACTER  STRING  WITHIN  THE 
PARENTHESIS  TO  AN  INTEGER  NUMBER  * 

IF <  JPAREN  ♦ EQ ♦  0  )  GO  TO  35 

IF (  IPAREN  ♦ NE ♦  2  )  THEN 

PRINT  %f'  ERROR  -  PARENTHESIS  MISSING  IN  THE  STRING  '  f  TEXT { £  ♦  I.  T ) 
PRINT  *, '  ' 

ISTAT  =  1 
GO  TO  9999 
ELSE 

CALL  ELEMNT (  TEXT ( 1 JLT )  *  L.T  f  IPOINT  f  ISTAT  ) 

IF <  ISTAT  »NE»  0  )  GO  TO  9999 
IPAREN  =  0 
ENPIF 

i  CONTINUE 


DETERMINE  WETHER  ANY  BLANKS  EMBEDDED  WITHIN  THE  QUANTITY ♦  IF  SO 
THEN  ITS  AN  ERROR ♦ 

N  =  INDEX <  TEXT (  KLT  )  f'  '  ) 

IF(  N  .;iE»  0  >  THEN 

PRINT  %t'  ERROR  -  CAN  NOT  HAVE  BLANKS  EMBEDDED  WITHIN  THE  NAME ' r 
%  TEXT  (  UI.T  > 

PRINT  :tr'  ' 

ISTAT  =  I 
GO  TO  9999 
END  IF 

SEARCH  FOR  THE  QUAHTJTV  IN  THE  LIST  NAMLIS ♦  IF  NOT  FOUND  -  ERROR 
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DO  100  I  =  1  iNVAR 

IF (  TEXT <  1JLT  )  *EQ.  NAMLIS<  J  )  )  00  TO  200 
100  CONTINUE 
C 

PRINT  Xt  '  ERROR  -  COULD  NOT  FIND  >  '  f  TEXT  (3  * LT> i '  <  IN  NAMELIST  ' 
PRINT  tt  '  ' 

I3TAT  =  1 
GO  TO  9999 
C 

200  CONTINUE 
C 

INL  =  I 
IEOUAL  =  0 
NUN  =  0 
IT EXT  =  0 
C 

IF  (  I  POINT  ♦  LE *  NDIM2<  INL  )  ♦  AND  *  IPOINT  » GE  •  NBIHX  (  INI  )  ) 

8  THEN 

1  TEXT  =  IPOINT  -  NDIM1 <  TNI.  ) 

IPOINT  =  1 
C 

ELSE 

C 

PRINT  *t'  ERROR  -  SPECIFIED  ARRAY  ELEMENT  OUT  OF  BOUNDS 
%  'FOR  THE  ARRAY  * ? NAMLTS (  INL  )  ♦  '  *' 

PRINT  X t'  THE  LOWER  AND  UPPER  BOUNDS  ARE  '> 

2  NDIM1C  INL  )  t  ;;DIM2<  INL  ) 

PRINT  Xt '  CHANGE  THE  ARRAY  ELEMENT  ' y IPOINT t '  TO  FIT'r 

S  '  WITHIN  THE  SPECIFIED  LIMITS*  ' 

ISTAT  =  1 
'GO  TO  999? 

C 

ENDIF 

FIND  THE  CHARACTER t INTEGER  OR  REAL  STRING  ANn  STORE  IN  THE 
APPROPRIATE  LOCATION  TO  BE  LATER  USED  BY  THE  MAIN  ROUTINE* 

210  CONTINUE 
C 

NJ  =  ISIZC  NDIMKJNL)  >  NDIM2CJNL)  ) 

C 

NTCI  =  NCHAR  +  NINTG 
C 

IF (  INL  » LE  *  NCHAR  )  THEN 
C 

POINTR  =  CPOSC  INL  ) 

CALL  FNDCHRC  STRING>ICHR(  POINTR  )  ?  I  TEXT  » L  S>  L  Tf  NL  r  NJ  r  IF  GOAL  t 
X  NUM r ISTAT  ) 

C 

IF (  ISTAT  *  EQ •  1  )  GO  TO  9999 

w 

ELSE  IF (  INL  ♦ GT  *  NCHAR  .AND.  INL  *i  F*  NTCI  )  THEN 
C 

II  =  INL  -  NCHAR 
POINTR  =  IPOS<  IT  ) 

CALL  END I NT (  STRING  t  INTGR(  POINTR  )  t  CDATA  *  I DAT  f  NDAT  t 

3  NIDAT  t  I  TEXT  t  LS  ♦  LT  ♦  NL  >  NJ  t  IEOUAL  t 

%  NUM  f  ISTAT  ) 

C 

IF (  ISTAT  *  EQ ♦  1  )  GO  TO  9999 
C 
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ELSE  IF  C  XNL.  *GT .  NTCI  )  THEM 
C 

II  -  INL  -  NTCI 
POINTR  =  RPQS (  II  ) 

CALL  FNDRL. (  STRING  *  REALS v  FGINTR  )  *  CDATA  *  RDAT  >  NDAT  , 

%  NR  OAT  *  ITEXT  *  1.3  *  LT  >  Ml.  *  NJ  *  IEQUAI.  * 

3  NUM  *  ISTAT  ) 

L 

IF  <  J3TAT  ♦ EG  *  1  )  GO  TO  9999 
ENDIF 

IF  THE  POINTER  IS  AT  THE  END  GF  THE  INPUT  STRING  x  ♦  * t  THE  PARSING 
OF  THE  INPUT  STRING  HAS  SEEN  COMPLETED  THEN  GO  TO  ????  TO  THE  MAIN 
ROUTINE*  OTHERWISE*  GO  TO  999?  TO  THE  BEGINNING  AND  SEARCH  FOR  THE 
NEXT  QUANTITY  *  ETC  * 

DO  215  I  =  1 f NTEXT 
215  TEXT  <  I : I  )  =  '  ' 

IF  (  NL  ♦  L  T *  LS  )  GO  TO  10 


IF  THE  SECOND  DOLLAR  SIGN  WAS  NOT  FOUND  THEN  INCREMENT  THE  LINE 
COUNTER <  •*..•>♦  OTHERWISE*  RESET  THE  LINE  COUNTER  TO  ONEQ)  AND  THE 
ECHO  PARAMETER  TO  ZEROCOK 


IF  <  NDS  ♦  Ffl  ♦  0  )  THEN 
L  =  L  +  1 
RETURN 
ENDIF 
C 

9999  CONTINUE 
C 

L  =  1 
IECHO  =  0 
ITEXT  a  0 
C 

RETURN 

END 


SUBROUTINE  FNDCHR!  STRING  >  X C  t  XTEXT  >  LS  >  L.T  r  L  *  NJ  > 

X  ,T EQUAL  >  NU»1  t  XSTAT  ) 

C 

C  THIS  SUBROUTINE  FINDS  A  CHARACTER  VALUE  AND  STORES  IT  JN  THE 
C  APPROPRIATE  ARRAY (  •IC1  ). 

C 

c  input: 
c 

C  STRING  - 

C 

C  LS  - 

C 

C  LT  - 

C 

c 

C  IC  - 

c 
c 

C  L  - 

C 
C 

C  NJ  - 

C 

C 

CHARACTER*!*)  STRING  >  IC 
CHARACTER*20  TEXT 

c 

PARAMETER!  NTEXT  »  20  ) 

C 

DIMENSION  ICC  NJ  ) 

C 

C**************** ********** ******** ****************** ******************* 
C*************** ******************************************************** 

c 

ICO  =  ICHAR (  'O'  ) 

IC?  =  ICHAR  <  ' 9 '  ) 

C 

DO  5  I  =  1 t NTEXT 
5  ‘  TEXT !  HI  )  «  '  ' 

C 

102  =  0 
C 

IL  =  L  -  LT  -  1 
C 

10  CONTINUE 
C 

DO  100  I  =  L  f LS 
C 

C  CHECK  FOR  BLANKS.  IF  SO  -  INCREMENT. 

C 

IF <  STRING <  X : J  )  .EG.  '  '  >  GO  TO  100 
C 

C  LOOK  FOR  AN  EQUAL  SIGN  AFTER  THE  QUANTITY .  IF  NOT  THERE  tiR  0 i HER 
C  CHARACTER  IS  THERE!  gXCUJQXHI*  BLANKS  >  THEN  XT'S  AN  ERROR. 

C 

IF!  STRING  i  X:i  )  .EG.  '  =  "  )  THEN 
I EG HAL  -1 
GO  TO  100 
EKOIr 


CHARACTER  STRING  OF  LENGTH  *  LS  * ♦ 

LENGTH  OF  CHARACTFR  STRING. 

LENGTH  OF  THE  QUANTITY  WHICH  IS  F3UIVALENCED  TO  THE 
ARRAY  'IC*. 

CHARACTER  ARRAY  CONTAINING  THE  INPUT  CHARACTER  VALUE. 
DIMENSIONED  TO  'NJ*. 

POINTER  IN  THE  STRING  ‘STRING'.  AT  THE  END  OF  THE 
ROUTINE  IT'S  UPDATED  TO  RETURN  THE  NEW  POSITION. 

DIMENSION  OF  THE  ARRAY  IC 


oooo  noooo  onooooo  o  on  o  n  o  o  n  o  n 


IF  EQUAL  SIGN  FOUND  THEN  NOW  SEARCH  AND  FIND  THE  CHARACTER  NEEDED  FOR 
THE  INPUT  * 

IF <  X EQUAL  ♦ Eft «  i  )  THEN 

IF  A  COMMA  FOLLOWS  AN  EQUAL  SIGN  THEN  IT'S  AN  ERROR ♦ 

IF  (  STRING  C  I : X  )  *EQ*  )  THEN 

IF(  ITEXT  *EQ  *  0  )  THEN 

PRINT  X?'  ERROR  -  COMMA  IN  WRONG  POSITION 
Z  STRING  <  ILM  )>' 

PRINT  Xt '  ' 

ISTAT  =  1 
RETURN 
END  IE 
GO  TO  100 
ENDIF 

IF (  STRING (  III  )  ♦ EQ ♦  )  THEN 

IF (  ITEXT  ♦ En »  0  )  THEN 

PRINT  *?'  ERROR  -  DOLLAR  SIGN  IN  WRONG  POSITION  '> 
l  STRING  (  Il.:X  )r'<-' 

PRINT  *> '  ' 

ISTAT  =  1 
RETURN 
'  ENDIF 

L  =  I 
RETURN 

ENDIF 

CHECK  FOR  A  QUOTE  -  INDJCATFS  CHARACTER  VALUE  AHEAD  AND  FIND  THE 
CHARACTER  STRING  USING  ROUTINE  *CHRVAL#* 

ICHS  =  ICHAR (  STRING (  151  )  ) 

IF  C  STRING  <  I J  J  )  ♦  EO  ♦  ""  )  THEN 

CALL  CHRVALI  STRING  tit  LS  >  IQ?  t  ISTAT  ) 

IF (  ISTAT  *  EQ  *  1  )  RETURN 

ML  =  JG2  -  J  -  1 

NUM  =  0  INDICATES  THAT  THE  CHARACTER  VALUE  WILL  NOT  Bn  MULTIPLIED  BY 
AN  INTEGER ♦  THERFORE*  CHARACTER  VALUES  WILL  BE  EXPECTED* 

IF  <  NUM  *  EQ  *  0  )  THEN 
ITEXT  =  ITEXT  +  1 

JTEXT  COUNTS  THE  NUMBER  OE  CHARCA7FR  VALUES*  IF  ITS  GREATER  THAN  THE 
SPECIFIED  DIMENSION  OF  THE  ARRAY  THEN  IT'S  AN  ERROR* 

IF  C  JTFXT  *  GT  *  NJ  )  THEN 

PRINT  *>'  ERROR  -  TOO  MANY  VALUES  FOR  '  t STRING <  XL *L  > 

PR IN7  ' 

ISTAT  =  1 
RETURN 
ENDIF 
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O  U  U 


ICC  ITFXT  )  =  STRING  C  > 

C 

r 

DC  20  J  =  ITXTO  >  ITEXT 

IC\  J  1  ~  STRING <  IM:iQ2-t  > 

20  CONTINUE 

C 

NUM  =  0 
ENDTF 


SET  NEW  POINTER  POSITION  WITHIN  "STRING  *  AND  GO  RACK  TO  THE  Fh GINNING 


L  =  IOC  f  1 
00  TO  10 


C  DETERMINE  Jr  A  CHARACTER  STRING  IS  PRECEDED  RY  AN  INTEGER  A NO  A 
C  MULTIPLICATION  SIGN  INDICATING  THAT  THE  CHARACTER  VALUE  WILL  HE 

C  STORED  IN  THE  ARRAY  *10"  THE  AMOUNT  GIVEN  SY  THE  INTEGER ♦ 

C 


C 


c 

c 


c 


c 


c 


c 

50 

C 


c 


ELSE  XF<  ICHS  ♦  GE  ♦  ICO  ♦AND*  ICHS  ♦!.£,  1C9  )  THEN 

IF  (  TEXT  (  .t  :.t  )  •  EQ  ♦  '  '  )  THEN 
TEXT <  HI  )  =  STRING <  XJX  ) 

JL  -  1 
ELSE 

TEXT  =  TEXT (  JJJL  )//STRING<  lit  ) 

JL  »  JL  +  1 
ENDIF 

ELSE  Jr(  STRING<  I J J  >  ♦ EC ♦  )  THEN 

IF <  TEXT C  1*1  )  « EH ♦  '  '  )  THEN 

PRINT  ERROR  -  NEED  AN  INTEGER  TO  MULTIPLY  'p 

*  CHARACTER  VALUE  '*STRING<  XI.  :Xr7  <~ ' 

PRINT  ' 

ISTAT  ^  1 
RETURN 
ELSE 

CALL.  CONVER  <  TEXT  (  1 J  JL  )  t  VAR  t  ISTAT  ) 

IP (  ISTAT  ♦ EU ♦  X  )  RETURN 

ITXTC  =  ITEXT  +  1 

ITEXT  =  ITEXT  +  INT(  VAR  ) 

IF (  ITEXT  ♦ 6T •  NJ  )  T HEN 

PRINT  *>'  ERROR-  TOO  MANY  VALUES  FOR  't STRING < IL t XV/7 
PRINT  ' 

ISTAT  =  1 
RETURN 
ENDIF 

NUM  =  1 

DO  50  J  -  .UNTEXT 
TEXT (  JJ.J  )  =  '  ' 


ENDIF 

ELSE  IF  (  TEXm:X)  *NE* 


*  AND*  STRING  <  IM  >  >NE,  )  THEN 
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PRINT  #>'  WARNING  -  POSSIBLE  ERROR  '>  STRING  v  II.  JX  >?'<-' 
PRINT  *j'  ' 

C 

L  =  I 
RETURN 

ELSE 

IF<  ,1  TP  XT  ♦  L  E  *  NJ  )  THEN 
L  =  I 
RETURN 
ENDIF 
C 

ENfilr 

C 

ELSE 

C 

PRINT  **'  ERROR  -  WRONG  CHARACTER  OR  EQUAL  SIGN  HISSING  '? 

3  STRING  C  XLtIFX  ) 

PRINT  #> '  ' 

ISTAT  *  1 
RETURN 
C 

ENDIF 

C 

100  CONTINUE 
r 

L  af  LS 

C 

RETURN 

END 
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SUBROUTINE  FNDINT (  STRING  »  IC  f  CDATA  ,  IDAT  r  NDAT  t  NIDAT  > 

%  TTSXT  .*  L3  »  LT  >  l.  .  NJ  t  IEQUAI.  f  NlJfi  , 

8  ISTAT  ) 

CXXXXmXMtXtXXWXXXXXMttM'MXXXMMM’WWMMXTWMMXWWW-XMMttXX**'** 

c 

C  THIS  SUBROUTINE  FINDS  AN  INTEGER  VALUE  FROM  THE  INPUT  CHARACTER 
C  STRING  ‘STRING*  *  THIS  ROUTINE  ALLOWS  THE  USER  TO  SPECIFY 

C  ALPHANUMERIC  CHARACTERS  FOR  THE  QUANTITY  WHICH  WILL  THEN  BE 

C  CONVERTED  TO  UNIQUE  INTEGER  VALUES. THIS  PART  IS  ACCOMPLISHED  IN 

C  SUBROUTINE  CCI. 

C 

c  input/output: 
c 

C  STRING  -  CHARACTER  VARIABLE  Or  I.ENTH  M.S‘  CONTAINING  THE 

C  INPUT  CHARACTER  STRING. 

c 

C  LS  -  LENGTH  GF  CHARACTER  STRING ♦ 

C 

C  LTEXT  -  LENT6HT  OF  THE  QUANTITY  WHICH  IS  EQUIVALENT  TO  THE 

C  ARRAY  -IC** 

C 

C  IC  -  INTEGER  ARRAY  CONTAINING  THE  APPROPRIATE  VALUES  FOR 

C  THE  SPECIFIED  QUANTITY* 

C 

C  L  -  POINTER  IN  THE  STRING  ‘STRING**  AT  THE  END  OF  THE 

C  ROUTINE  ,l*  IS  UPDATED  TO  RETURN  THE  NEw  POStTION* 

£ 

C  NJ  -  DIMENSION  OF  THE  ARRAY  ‘IC4* 

C 

C  CDATA  -  CHARACTER  ARRAY  CONTAINING  THE  CHARACTER 

C  REPRESENTATION  OF  VARIOUS  VALUES* 

C  CDATA  DIMENSIONED  TO  • NBA7* ♦ 

C 

C  IDA  r  -  INTEGER  ARRAY  CONTAINING  THE  RESPECTIVE  XfyT&G&s? 

C  VALUES*  DIMENSIONED  TO  ‘NIOAT", 

C 

C  NDAT  -  TOTAL  NUMBER  OE  VALUES <  INTEGER  AND  REAL  )  IN  THE 

C  ARRAY  -CBATA*. 

C 

C  NIDAT  -  TOTAL  NUMBER  OF  INTEGER  VALUES  IN  THE  ARRAY  "I DAT 1  * 

C 

C 

CHARACTERS*)  STRING  t  CDATA 
CHARACTER*?.*)  TEXT 
C 

PARAMETER  <  NTEX  f  =  -20  ) 

C 

DIMENSION  ICC  NJ  )  ,  CDATAC  NDAT  >  >  IDATC  NIDAT  > 

C 

C^******************************^*******************^***^*************** 

C********************ic**4:4c**.,t**5*;***********‘i:*******x.ic**)ky*4t^'«i(^i**^******. 

C 

CO  5  I  JjNTEXT 
S  TEXT <  til  >  S  '  ' 


IC.O  ~  ICHAR(  'O'  ) 
IC?  =  ICKAP.f  ) 

IGMS  =  !CHAf!<  ) 

ICPS  *  ICUARf  '  +  '  ) 
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IQ2  =  0 
C 

IL  =  L  -  LT  -1 
C 

10  CONTINUE 
C 

DO  100  I  =  LfL3 
C 


IF  < 

STRING  < 

x :  x  > 

.EG. 

'  '  )  GO  TO  300 

IF  ( 

STRING  < 

nz  ) 

♦  EG. 

>  THFN 

IEGUAL  =  1 
GO  TO  100 
ENDIF 


IF  THE  EQUAL  SIGN  HAS  BEEN  FOUND  THEN  FIND  THF  INTEGER  VALUES « 

IF <  IEGUAL  .  EG  *  1  )  THFN 

DETERMINE  IF  THE  CHARACTER  IS  A  NUMBER  *  IF  SO  y  C0NCA70NA7E. 

ICHS  =  JCHAR <  STRING <  TJX  )  ) 

IF (  {  ICHS  . GE ♦  ICO  *  AND  *  ICHS  * LF ♦  IC?  )  .OR* 

X  ICHS  » EG »  ICPS  .OR.  ICHS  .EG .  ICMS  )  THEN 

IF <  TEXT (  3 : 1  )  .EG*  '  '  )  THEN 
TF*'T (  i:i  )  =  STRING<  Itl  ) 

J  •  1 
ELSE 

TEXT  =  TEXT (  1JJL  )//STRTNGC  I J  X  ) 

JL  =  Jl.  +  1 
ENDIF 

CHECK  FOR  COMMA*  IF  NOTHING  IN  TEXT  -  ERROR .  OTHFRWISE  t  CONVERT  STRING 
TO  INTEGER. 


C 

C 

C 

C 

C 

s 

X 


ELSF  IF (  STRING (  III  )  .EG.  .OR.  STRING (  III  )  .EG.  ) 

THEN 

IF <  JG2  .EG.  9999  )  THEN 
IG2  =  0 

IF (  STRING (  J;i  )  .FG.  't '  )  GO  TO  100 

IF (  STRING (  XJI  )  .EG.  )  THEN 

L  =  I 
RETURN 
ENDIF 

ENDIF 

IF  (  TEXT  <  311  )  .EG.  '  '  .AND.  STRING  (  IJX  )  »nG.  '»'  ) 

THEN 

FRJNT  #>'  ERROR  -  NO  VALUE  FOUND  PRIOR  TO  REACHING  A  ' > 
'COMMA  't STRING (  ILJI  ) r ' 

PRINT  %$’  ' 

ISTAT  =  1 
RETURN 
ENDIF 
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IF  <  TEXT  (  an  )  .ECK  '  '  ♦  AND  ♦  STRING  (  IJI  )  ♦  EG  ♦  'I'  ) 

X  THEN 

PRINT  *>'  ERROR  -  NO  VALUE  FOUND  PRIOR  TO  REACHING  A  't 
S  'DOLLAR  SION  'f STRING <  ILJI  >;'  <-' 

PRINT  4C r  '  ' 

ISTAT  =  1 
RETURN 
END  IF 

NUN  =  0  INDICATES  NO  MULTIPLICATION  OF  VALUE  BY  INTEGER ♦ 

IF  <  NUN  ♦ FQ .  0  )  THEN 


I  TEXT  -  I  TEXT  *r  1 
1F<  I  TEXT  « GT »  N  J  )  THEN 

PRINT  *>'  ERROR  -  TOO  MANY  VALUES  FOR  'rSTRINOC  ILtl  > 
PRINT  %7  '  # 

ISTAT  =  1 
RETURN 
ENDIF 

CALL  CONVERC  TEXT <  UJL  )  >  VAR  t  ISTAT  ) 

IF<  ISTAT  » EG •  1  )  RETURN 

ICC  I TEXT  )  a  I NT (  VAR  ) 

DO  20  J  -  ItNTEXr 
20  TEXT  <  UM  )  »  '  ' 

ELSE 

CALL  CONVERC  7 EXT (  J Ml  >  >  VAR  ?  ISTAT  ) 

IF <  1ST AT  ♦ EG ♦  1  )  RETURN 

I VAR  =  I NT <  VAR  ) 

DO  30  J  =  JTXTO t TTEXT 
ICC  J  )  .*  I  VAR 
30  CONTINUE 

NUM  *  0 


DO  40  J  *  1 t NTEXT 

40  TEXT <  JtJ  )  ~  ' 

END*  F 


RESET  POINTER  AND  GO  TO  BEGINNING  TO  FIND  ANOTHER  VALUE \  \r  ANY  ) 

IF<  STRING  <  IJI  )  *E0*  5  THEN 

l  =  r 

RETURN 

ENDIF 

L  p  I  +  1 
GO  TO  10 

IF  **a  FOUND  THEN  CONVERT  TEXT  TO  INTEGER  AND  MULTIPLY  BY  THE  UPCOMING 
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VAI.IJE » 

ELSE  IF <  STRING?  X!I  )  .EQ.  >  THEN 

IF  (  TEX  f  ?  XU  )  ,EQ.  '  '  )  THEN 

PRTNT  *>'  ERROR  -  INTEGER  REOIIIRE!)  '  >  STRING?  Xl.U  )> 

1 

ISTAT  =  1 
RETURN 
END  IF 

KL  -  JLEN<  TEXT  >  N  TEXT  ) 

N  INDEX <  TEXT <  1IKL  >  >  > 

i F  <  N  » NF *  0  >  THEN 

PRINT  '  ERROR  -  CAN  NOT  MULTIPLY  BY  NEGATIVE  INTEGER  '? 
;i  TEXT  (  tJKL  ) 

PRINT  *F'  ' 

ISTAT  =  1 
RETURN 
ENDIF 

CALL  CONVER <  TEXT <  UKl  )  r  VAR  t  ISTAT  ) 

IF C  ISTAT  .EG#  1  )  RETURN  - 

ITXTO  =  ITEXT  +  1 

ITEXT  =  ITEXT  +  INT<  VAR  ) 

IF (  ITEXT  * GT ♦  NJ  )  THEN 

PRINT  *>'  ERROR  -  TOO  MANY  VALUES  FOR  '> STRING <  XL S X  ) 
PRINT  *7 '  ' 

ISTAT  =  1 

Return 

enbif 

NUM  =  1 

DO  *0  J  =  IjNTEXT 
60  TEXT <  JJJ  )  *  '  ' 

IF  QUOTE  FOUND  THEN  FIND  THE  CHARACTER  VALUE  CONVERT  IT  TO  A  UNIQUE 
NUMBER  USING  ROUTINE  *CCI-> 


C 

C 

C 


C 


C 

c 

c 


ELSE  IF  (  STRING  (  j:j  )  ♦  EQ  ♦  ""  )  THEN 

CALI.  CHRVAL<  STRING  i  I  t  LS  t-  102  y  ISTAT  ) 

IF (  ISTAT  ♦ EO ♦  1  )  RETURN 

CALI.  CCI  (  STRING  (  JJIQ2  )  j  IVAR  y  CD  AT  A  y  IDAT  t  NDAT  y 

NIDAT  f  ISTAT  ) 

IF (  ISTAT  *En*  X  )  RETURN 
IF (  HUH  • EQ  >  0  )  THEN 
I  TEX  F  ^  ITEXT  +  1 
IF (  ITEXT  *6T •  NJ  )  THEN 

PRINT  *?'  ERROR  -  TOO  MANY  VALUES  FUR  * y STRING!  sllS  > 
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PRINT  ' 

ISTAT  =  1 
RETURN 
ENDIK 

IC<  I TEXT  )  -  IVAR 


C 


30 

C 


ELSE 

DO  80  J  =  I TXTO  ? I TEXT 
IC(  J  )  =  I VAR 
CONTINUE 

HUM  =  0 


C 

C 

c 


END  IF* 

L  =  IQ?  +  i 
IQ?  =  9??? 

GO  TO  30 

ELSE  IF <  TEXT ( X  •  X )  ♦ NF .  '  '  * AND  *  STRING (  III  )  *NE «  )  THEN 


C 

C 

c 

c 

c 

c 

c 

100 

c 

c 


PRINT  *>'  WARNING  -  POSSIBLE  ERROR  'jSTRINGC  XI.  1 1  >*'<-' 
PRINT  *7/  ' 

L  *  I 
RETURN 

ELSE 

IF <  ITEXT  J.E»  NJ  )  THEN 
L  a  I 
RETURN 
ENOIF 

Cv* 

ENDIr 


ELSE 

PRINT  ERROR  -  WRONG  CHARACTER  OR  EQUAL  SIGN  MISSING  '7 

Z  STRING <  ILM  )$’  <-' 

PRINT  *7 '  ' 

ISTAT  =  1 
RETURN 

END  IF 
CONTINUE 


RETURN 

END 
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SUBROUTINE  FNDRL<  STRING  >  IC  >  CDATA  >  RDAT  >  NDAT  f  NRDAT  ? 
t  ITEXT  f  LS  t  LT  r  I.  t  NJ  f  I  EQUAL  t  NUN  * 

&  ISTAT  ) 

C 

C  THIS  SUBROUTINE  FINDS  AN  INTEGER  VALUE  FROM  THE  INPUT  CHARACTER 
C  STRING  *  STRING “ *  THIS  ROUTINE  ALLOWS  THE  USER  fO  SPECIFY 

C  ALPHANUMERIC  CHARACTERS  FOR  THE  QUANTITY  WHICH  WIU  THEN  BE 

C  CONVERTED  TO  UNIQUE  INTEGER  VALUES*THI3  PART  IS  ACCOMPLISHED  IN 

C  SUBROUTINE  CCI. 

C 

c  input/output: 

c 

C  STRING  - 

c 
c 

C  LS  - 

C 

C  LT EXT  - 

C 
C 

C  IC  - 

C 
C 

C  L  - 

c 

c 

C  NJ  - 

C 

C  CDATA  - 

C 
C 
C 

C  RDAT  - 

C 

c 

C  NDAT  - 

C 
C 

C  NRHAT  - 

C 

C#*#****#***###*##**#***#*#*###*#**#**##***#*#*****#*#********#****#*#** 

C 

CHARACTER# < t )  STRING  >  CDATA 
CHARACTER*?*)  TEXT 
C 

REAL  IC 
C 

PARAMETER <  NTFXT  =  20  ) 

C 

DIMENSION  ICC  NJ  )  >  CDATA<  NDAT  )  >  RDAT <  NRDAT  ) 

C 

C**  *************************************<**********  ****,i';f::i'*****#********* 

c 

DO  5  I  =  U NTFXT 
5  TEXT (  t:t  /  *  '  ' 

C 

ICO  =  JCHAR(  'O'  ) 

IC?  =  I CHAR (  '9'  ) 

I CMS  =  ICHAR<  '-'  ) 


CHARACTER  VARIABLE  OF  LENTH  *  LS  *  CONTAINING  THE 
INPUT  CHARACTER  STRING. 

LENGTH  OF  CHARACTER  STRING ♦ 

LENTGHT  OF  THE  QUANTITY  WHICH  IS  EQUIVALEND  TO  THE 
ARRAY  VtC. 

INTEGER  ARRAY  CONTAINING  THE  APPROPRIATE  VALUES  FOR 
THE  SPECIFIED  QUANTITY  * 


POINTER  IN  THE  SIRING  "STRING * ♦  AT  THE  END  OF  THE 
ROUTINE  *  L #  IS  UPDATED  TO  RETURN  THE  NEW  POSITION* 

DIMENSION  OF  THE  ARRAY  "  JC"  * 

CHARACTER  ARRAY  CONTAINING  THE  CHARACTER 
REPRESENTATION  OF  VARIOUS  VALUES. 

CDATA  DIMENSIONED  TO  "NDAT* ♦ 

REAL  ARRAY  CONTAINING  THE  RESPECTIVE  REAL 
VALUES.  DIMENSIONED  TO  "NRDAT • ♦ 

TOTAL  NUMBER  OF  VALUES  (  INTEGER  AND  REAL.  )  IN  THE 
ARRAY  "CDATA'. 

TOTAL  NUMBER  OF  REAL  VALUES  TN  THE  ARRAY  ' RDAI " ♦ 


88 


oooo  o  oonoooon 


ICPS  «  ICHAR (  '  +  '  ) 
ICPER  ^  I  CHAR  C  ' * '  ) 


IQ  2  =  0 
C 

IL  a  L  -  LT  -  1 
C 

10  CONTINUE 

i; 

DO  100  I  =  I.  f  I..S 

IF (  STRING*  JJI  >  .EQ.  '  '  >  GO  TO  100 

IF <  STRING  <  I : I  )  » EQ  *  '  =  '  >  THFN 
I EQUAL  =  1 
GO  TO  100 
ENG  IF 

IF  THE  EQUAL  SIGN  HAS  BEEN  FOUND  THFN  FIND  THE  REAL  VALUES  * 

IF <  JEQUAL.  *  FQ  ♦  1  )  THEN 

DETERMINE  IF  THE  CHARACTER  IS  A  NUMBER  *  IF  SOt  DETERMINE  WHETHER  IT'S 
A  ZERO*  IF  SO »  ERRORJOTHERWISE  CONCATONATE. 

ICHS  *  ICHAR (  STRING (  J!I  )  ) 


IF <  (  ICHS  . GE . 

ICO  .AND. 

ICHS  .LE. 

IC9  )  .OR. 

£ 

ICHS  .EQ. 

ICHS  *  OR  * 

ICHS  .EQ. 

ICPS  .OR. 

l 

ICHS  .EQ* 

ICPER  )  THEN 

IF <  TEXT*  X J 1 

)  .EQ.  ' 

'  )  THEN 

TEXT*  IJI  ) 

a  STRING* 

IJI  ) 

JL  =  1 
ELSE 

TEXT  =  TEXT (  3 l JL  )//STRING*  IJI  ) 

JL  a  JL.  +  1 
END  IF 

CHECK  FOR  COMMA.  IF  NOTHING  IN  TEXT  -  ERROR ♦ OTHERWISE  t CONVERT  STRING 
TO  REAL. 

ELSE  IF (  STRING <  IJI  )  . EO .  '?'  .OR.  STRING (  IJI  )  .EH.  ) 

3  THEN 

C 

IF <  102  *EQ .  9999  )  THEN 
C 

102  =  0 

c 

IF (  STRING <  IJI  )  .EG.  )  GO  TO  1O0 

IF  <  STRING  <  I : I  )  .EQ.  )  THEN 

L  *  I 
RETURN 
ENDIF 
C 

ENDIF 

C 

IF  (  TEXT  <  3U  )  *  Ffl .  '  '  .AND.  STRING  <  IJI  )  *  EQ  *  ) 

*  THEN 

PRINT  *t'  ERROR  -  Na  VALUE  FOUND  PRIOR  TO  REACHING  A  '> 
l  'COMMA  STRING*  XI.  S I  >>'  <-' 
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ISTAT  =  i 
RETURN 
ENOIF 
C 

IF  C  TEXT  <  U1  )  ♦  EG  ♦  '  '  » AND *  STRINGC  X :  J  )  ♦EH.  ) 

2  THEN 
PRINT  *»'  ERROR  -  NO  VALUE  FOUND  PRIOR  TO  REACHING  ft  'y 

3  'DOLLAR  SION  *,  STRING  <  XU  I  >;' 

PRINT  Xt '  ' 

ISTAT  =  1 
RETURN 

ENHIF 

NUN  =  0  INDICATES  NO  MULTIPLICATION  OF  VALUE  BY  XNTSGfc  E ♦ 

i'F  <  NUN  *  EQ  ♦  0  )  THEN 

ITEXT  =  XTEXT  +  1 
C 

IF  <  XTEXT  ♦OT,  N  J  )  THEN 

PRINT  *y '  ERROR  -  TOO  MANY  VALUES  FOR  '  y  STRING  <  XLJX  > 
PRINT  *y'  ' 

ISTAT  3  1 
RETURN 
ENDIF 
C 

CAIJ.  CONVERC  TEXT <  .t  i  JL  )  y  VAR  y  ISTAT  ) 

C 

IF <  ISTAT  ♦ EG ♦  1  )  RETURN 
IC<  XTEXT  )  =  VAR 
DO  00  J  *  3  yNTEXT 

20  TEXT  C  JLJ  )  =  '  ' 

C 

ELSE 


C 

C 


SO 

c 

c 

40 


CALI.  CONVERC  TEXT (  X  *  JL  )  y  VAR  y  ISTAT  ) 

IF (  ISTAT  ♦ EG  *  3  )  RETURN 

DO  30  J  3  XTXTO  y  X  TEXT 
ICC  J  )  =  VAR 
CONTINUE 

NUN  =  0 

DO  40  J  -  1 y NTE XT 
TEXT C  JJJ  )  =  '  ' 


ENDTF 

C 

C  RESET  POINTER  AND  00  TO  BEGINNING  TO  FIND  ANOTHER  VALUE <  IF  ANY  > 
C 

IFC  STRING  C  X  J I  >  .EG.  'S'  )  THEN 
L  -  I 
RETURN 
E«niF 
c 

L  =  X  +  1 
GO  TO  .10 
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non  o  o  noon  no  n  o  o  non 


UPUO MING 


*  fc '  FOUND 
VALUE  » 


THEN  CONVERT  TEXT  TO  INTEGER  AND  MULTIPLY  J5Y  THE 


ELSE  Jr  <  STRING <  III  >  ♦EO.  'S'  )  THEN 
IF-;  TEXT (  1 : 1  )  *EQ*  '  '  )  THEN 

PRINT  %j'  ERROR  -  INTEGER  REQUIRED 

PRINT  *t  '  ' 

ISTAT  “  1 
RETURN 
END  IF 


KL  =  X !.  F N <  TEXT  ;  NTEXT  > 

N  =  INDEX (  TEXT \  J SKL  )  ?  ) 


STRING (  XL  T  X 


IF<  N  *  NE ♦  0  )  THEN 

PRINT  *>'  ERROR  -  CAN  NOT  MULTIPLY  SY  NEGATIVE  INTEGER  '> 
TEXT <  i:kl  ) 

PRINT  *r'  7 
ISTAT  -s  1 
RETURN 
ENDIF 


N  *  INDEX (  TEXT <  JUKI.  )>'♦') 


IF <  N  ♦ NE ♦  0  )  THEN 

PRINT  *>'  ERROR  -  LOOKING  FOR  AN  INTEGER  NOT  A  REAL  '> 
X  'NUMBER  't  TEXT <  JUKI.  ) 

PRINT  *>'  ' 

ISTAT  *  1 
RETURN 
ENDIF 

CALL  CONVER (  TEXT (  XJKL  )  t  VAR  t  ISTAT  ) 

IF <  ISTAT  ♦ EO ♦  X  )  RETURN 

ITXTO  =  IT EXT  +  l 

ITEXT  =  IT  EXT  +  INK  VAR  ) 

IF (  ITEXT  *GT 4  NJ  )  THEN 

PRINT  *t'  ERROR  -  TOO  MANY  VALUES  FOR  ' > STRING!  ILJX  ) 
PRINT  * 7 '  ' 

ISTAT  =  1 
RETURN 
ENDIF 

NUrt  =  1 


DO  60  J  =  tf NTEXT 
60  TEXT <  JJJ  )  =  '  ' 

IF  QUOTF  FOUND  THEN  FIND  THE  CHARACTER  VALUE  CONVERT  IT  TO  A  UNIQUE 
NUMBER  USING  ROUTINE  #CC Rg, 

ELSE  IF <  STRING <  j;i  )  ♦EO*  ""  )  THEM 

CALL  CHRVAL '  STRING  r  I  r  LS  >  102  >  ISTAT  ) 


C 

c 

c 
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IF <  1ST AT  .FQ.  1  )  RETURN 

CALL  OCR  <  STRINGC  XSX&3  )  ?  VAR  t  CBATA  ? RDnT  r  ? 

NROAT  f  ISTAT  ) 

IF  C  I ST AT  .EG.  3  >  RETURN 
IF (  NUN  »EQ.  0  >  THEN 
ITEXT  =  IT EXT  r  1 
Iri  ITEXT  .GT.  NJ  )  THEN 

PRINT  Mr'  ERROR  -  TOO  MANY  VALUES  FOR  '.-STRING!  XI.  :  a  ) 
PRINT  *; '  ' 

ISTAT  =  i 
RETURN 
END  IF 

ICC  ITEXT  >  =  VAR 
ELSE 

DO  SO  J  =  ITXTO* ITEXT 
IC<  J  )  =  VAR 
CONTINUE 

NUN  =  0 

END  IF 

L  =•  IQ  2  +  1 
IQ?.  =  9??? 

60  TO  10 

ELSE  IF (  TEXTCIJI)  .NE.  '  '  .AND.  STRING!  JJX  >  ♦ NF ♦  'S'  )  TKh M 

PRINT  1 1  '  WARNING  -  POSSIBLE  ERROR  '  r STRING!  XI.  J  X  >?'<-' 
PRINT  *7'  ' 

L  *  I 
RETURN 

ELSE 

IF <  ITEXT  .LE.  NJ  )  THEN 
L  =  I 
RETURN 
END  IF 

ENDJF 

ELSE 

PRINT  Mr'  ERROR  -  WRONG  CHARACTER  OR  FOUAI.  SIGN  nXSSING  'r 

STRING*;  It.  .1  )*'  “C- ' 

PRINT  Mr'  ' 

ISTAT  =■  1 
RETURN 


END  IF 


100  CONTINUE 
C 

L  =  US 
C 

RETURN 

END 
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o  o  o  n  n  o  n  o  o  o  ri 


SUBROUTINE  CHRVAL (  STRING  ,  .1  >  LS  r  IQ?  r  ISTAT  ) 

CHARACTERS <  *)  STRING 

THIS  SUBROUTINE  FINDS  A  CHARACTER  VALUE  CHECKS  IT  AND  MAKES  SURE 
THAT  IT  IS  CORRECT »  IT  RETURNS  THE  VALUE  OF  THE  SECOND  QUOTE* 

MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM 

NC  -  LS  -  1 

IQ?  =  INDEX (  STRING (  I  +  IJNC  )  r  ""  ) 

IF <  IG2  ♦ EQ ♦  0  )  THEN 

PRINT  Mr  'ERROR  -  MISSING  SECOND  QUQTB  STRING <  UNO  )? 

3  '  <- ' 

ISTAT  =  1 
RETURN 
ENDIF 

CHECK  FOR  EMBEDDED  BLANKS  WITHIN  THE  TWO  QUOTES. 

IQ?  a  IQ?  ■*■  I 

N  a  INDEX<  STRINGC  i:iQ2  )  r  '  '  ) 

IF  <  N  *  NF.  *  0  )  THEN 

PRINT  Mr'  ERROR  -  CANNOT  HAVE  EMBEDDED  FLANKS  'j 
S  'IN  CHARACTER  VALUES  '> 

l  STRING (  IIIQ2  )  r'  <-' 

ISTAT  *  1 
RETURN 
ENDIF 
C 

RETURN 

END 
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SUBROUTINE  CONVERC  TEXT  t  VAR  t  ISTAT  ) 
C**XXXX*X%X*X%*X**X**X**X%*X**X******t**XX***t**Xt%*X****X***X**X**X*XXX 

c 

C  THIS  SUBROUTINE  WILL  CONVERT  AN  ANSI  REPRESENTATION  OF  A  NUMBER  TO 
C  AN  ACTUAL  NUMBER <  INTEGER  OR  REAL  -  EXCLUDING  1 E *  FORMAT  )* 

C 

C  INPUT/OUTPUT: 

C 

C  TEXT  -  CHARACTER  STRING  WITH  1 HF  ANSI  RFPRE3ENTATI0N  OF  A 

C  NUMBER 

C 

C  VAR  -  ACTUAL  NUMERICAL  VALUE  BEING  RETURNED ♦  IT  CAN  BE  EITHER 

C  INTEGER  OR  REAL  *  FOR  THE  MOMENT  A  REAL  NUMBER 

C  REPRESENTED  BY  *E*  FORMAT  IS  NOT  AVAILABLE* 

C 

C  >»»N0TE:  -  THE  ALGORITHM  FOR  THIS  ROUTINE  PROVIDED  BY 
C  FRANK  MANSFIELD  OF  N3WC*CA* 

C 

C 

CHARACTER #<*)  TEXT 
INTEGER  DIGIT 
C 

VAR  =  0*0 
C 

ICO  =  ICHAR (  '0'  ) 

C 

N  =  LEN <  TEXT  ) 

NT  =  ILEN <  TEXT  r ' N  ) 

C 

DO  10  I  =  1  >  NT 
C 

IF (  TEXT (  X : T  )  *  EQ ♦  '  '  )  GO  TO  10 
C 

IF <  TEXT (  I  *  I  )  *E0.  )  THEN 

C 

NSTR  =1+1 
SIG  =  -1*0 
C 

ELSE  IF <  TEXT <  ITT  )  .EQ*  '  +  '  )  THEN 
C 

NSTR  =1+1 
SIG  =  t.O 


C 

C 

C 

C 

C 


ELSE 

NSTR  =  I 
SIG  =  1.0 

ENDIF 

GO  TO  20 


10  CONTINUE 
C 

PRINT  ERROR  -  NO  CHARACTER  STRING  FOUND  JN  TEXT  TEXT  CUNT  ) 

ISTAT  =  1 
RETURN 
C 

20  CONTINUE 
C 
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NDEC  =  INDEX (  TEXT <  NSTRJNT  ) 

C 

IF<  NDEC  *EQ *  0  )  THEN 

NDEC  =  NT+1 

ELSE 

NDEC  =  NDEC  +  NSTE  -  1 

ENDIr 

C 

TEN  =0*1 
C 

DO  100  J  =  NDEC-1 rNSTR'-l 
C 

DIGIT  =  ICHhR (  TEXT (  III  )  ) 
TEN  =  JLO»*TEN 
VAR  =  VAR  +  TENKDIGIT 
C 

100  CONTINUE 
C 

TEN  =1,0 
C 

DO  200  I  =  NDEC+lfNT 
C 

DIGIT  =  ICHAR (  TEXT <  i:i  )  ) 
TEN  =  0*1*TEN 
VAR  =  VAR  +  TEN*DJGIT 
C 

200  CONTINUE 
C 

VAR  =  VAR*SIG 
C 

RETURN 

END 


t  '  ♦  '  ) 


-  ICO 


-  ICO 
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SUBROUTINE  CCI<  TEXT  >  JVAR  >  CDATA  t  IDAT  >  NBAT  ,  NIDAT  >  ISTAT) 
************)M** ************** ***************************************** 

THIS  SUBROUTINE  CONVERTS  CHARATER  INFORMATION  INTO  UNIQUE  INTEGER 
VALUES  BASE  UPON  THE  ALLOWABLE  INPUTS  FOR  THE  DESIRED  QUANTITY » 


C*********************************************************************** 

C 

CHARACTERS  <  *)  TEXT  >  CDATA 
C 

DIMENSION  CDATA  <  NDAT  )  >  IDAT <  NIDAT  ) 

C 

N  =  LEN (  TEXT  ) 

NT  =  I LEN C  TEXT  >  N  ) 

C 

DO  10  I  -  1 ; NBA T 

Xr<  TEXT <  i+i:NT-l  )  ♦ EQ ♦  CDATA <  I  )  >  GO  TO  20 
10  CONTINUE 
C 

PRINT  %f*  ERROR  -  ' t TEXT (  it NT  >>'  NOT  FOUND  IN  CDATA  ARRAY  ' 

ISTAT  =  1 

RETURN 

C 

20  CONTINUE 
C 

IF (  I  .GT.  NIDAT  )  THEN 

PRINT  *>'  ERROR  -  INCORRECT  CHARACTER  STRING  ,  EXPECTING  't 
S  ' ONE  OF  THE  >E  ' 

PRINT  *r<  CDATA (  J  )  rJ  =*.t>NIDAT  ) 

PRINT  *t'  INSTEAD  RECEIVED  THIS  *  f TEXT (  1 ♦ NT  ) 

ISTAT  =  1 
RETURN 

ENDIF 

C 

I VAR  =  IDAT (  X  ) 

C 

RETURN 

END 
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SUBROUTINE  CCR<  TEXT  »  VAR  >  CDATA  >  RDAT  ?  NDAT  >  NRDAT  t  1ST AT  ) 
C*******************#******#************************r***** ************** 
c 

C  THIS  SUBROUTINE  CONVERTS  CHARATER  INFORMATION  INTO  UNIQUE  REAL 
C  VALUES  BASED  UPON  THE  ALLOWABLE  TNPUTS  FOR  THE  DESIRED  QUANTITY* 

C 

C*********************************************************************** 

c 

CHARACTER* HO  TEXT  r  CDATA 
C 

DIMENSION  CDATA  <  NDAT  )  >  RDAT (  MRDAT  ) 

C 

N  =  LEN<  TEXT  ) 

NT  “  ILEN (  TEXT  >  N  ) 

r 

DO  j  0  I  -  X } NDAT 

IF  (  TEXT (  l-MJNT-j  )  ♦  FQ .  CDATA  (  I  )  )  60  TO  70 
10  CONTINUE 
C 

PRINT  *> '  ERROR  -  *t TFXT<  1 INT  )r'  NOT  FOUND  IN  CDATA  ARRAY  ' 

ISTAT  =  1 
RETURN 
C 

20  CONTINUE  ' 

C 

NIDAT  *  NDAT  -  NRDAT 
IF  C  I  *LE.  NIDAT  )  THEN 
C 

N  =  NIDAT  +  1 

IF (  NIDAT  *  EQ »  NDAT  )  N  *  NDAT 
C 

PRINT  #»'  ERROR  -  INCORREXT  CHARACTER  VALUE  *  EXPECTING  '  r 
&  'ONE  OF  THESE  ' 

PRINT  *t(  CDATA (  J  )  tJ  -N?  NDAT  ) 

PRINT  *?'  INSTEAD  RECEIVED  THIS  '>TFXT(  1JNT  ) 

ISTAT  *  i 
RETURN 
END  IF 
C 

II  »  I  -  NIDAT 
VAR  -  ROAT(  II  ) 

C 

RETURN 

END 
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SUBROUTINE  ELEMNT (  TEXT  »  L.T  ,  IVAR  »  ISTAT  ) 

C******jC*  ************************************************************** 

c 

C  THIS  ROUTINE  CONVERTS  A  CHARACTER  STRING  WITHIN  A  SET  OE  PARENTHESIS 
C  TO  AN  INTEGER  W 

C  ARRAY ♦ 

r* 

W 

c  infut/output: 

c 

C  TEXT  - 

C 
C 

C  LT  - 

C 

C  I VAR  - 

C 
C 
C 

C  ISTAT  - 

*C 

c 

CHARACTER*# ( * )  TEXT 
CHARACTER*.*,?  NUM 
C 

DO  1  I  =  1  >32 
1  NUM (  IJt  )  =  '  7 

C 

ICO  =  ICHARC  'O'  ) 

IC?  =  ICHAR (  '?'  ) 

C 

I CHIN  =  ICHAR  <  '-'  ) 

XCPLS  =  ICHAR  (  '  +  '  ) 

LOCATE  POSITION  OF  THE  SET  OF  PARENTHESIS*  THEN*  USE  THESE  VALUES 
IN  THE  DO  LOOP  TO  FIND  THE  CHARACTER  STRING  -  NUJi  *  WHICH  WILL  THEN 
BE  CONVERTED  TO  AN  INTEGER  VALUE  #IVAR*  ) 

NL  =»  INDEX (  TEXT  >  '('  ) 

NR  =  INDEX (  TEXT  t  ')'  ) 

C 

LT  =  NL  -  1 
C 

DC  10  I  a  Nl.  +  I  t  NR-1 
C 

IF <  TEXT (  IJT  )  *  EO ♦  '  '  )  GO  TO  10 
C 

ICHS  =  ICHAR (  TEXT (  J  J T  )  ) 

C 

IF (  <  ICHS  »G£ ♦  ICO  ♦  AND ♦  ICHS  .  I.E »  IC?  )  OR. 

S  ICHS  .EQ »  I CM IN  .OR.  ICHS  .EG.  XCPLS  )  THEN 

C 

IF  <  NIJH  <  1J1  )  *  EO »  '  '  )  THEN 
NUM (  111  )  =  TEXT (  X  J I  ) 

IN  =  1 
ELSE 

NUM  =  NUM<  III H  ) //TEXT (  111  ) 

IN  =  IN  +  1 
ENDIF 
C 

ELSE 


INDICATING  THE  POSITION  WITHIN  THE  SPECIFIED 


CHARTACTER  STRING  CONTAINING  THE  ARRAY  NAME  AND 
THE  ARRAY  ELEMENT. 

CHARACTER  LENGTH  OF  THE  VARIABLE. 

INTEGER  NUMBER  (  POSITIVE  OR  NEGATIVE  )  BEING 
RETURNED  TO  THE  CALLING  ROUTINE.  INDICATING  THE 
ARRAY  ELEMENT  DESIRED. 

ERROR  FLAG  ?  0  -  NO  ERROR?  1  -  ERROR 
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PRINT  *> 


ERROR  -  INCORRECT  CHARACTER  WITHIN  PARENTHESIS 
CHECK  'r  TEXT (  t  *  I  )t'  <=' 


PRINT  Xf  '  ' 

I3TAT  *  1 
RETURN 

END  IF 

CONTINUE 

CALL.  CONVF.R  (  KCJH  >  VAR  >  ISTAT  > 
IF <  ISTAT  ♦ NE ♦  0  )  RETURN 
I  VAR  =  INK  VAR  ) 


RETURN 

END 


FUNCTION  ILEN<  TEXT  t  N  ) 

CHARACTER* ( * )  TEXT 
C 

00  100  I  = 

IF  <  TEXT <  It  I  )  #  NE ♦  '  '  )  GO  TO  200 
100  CONTINUE 
C 

ILEN  0 

PRINT  *>'  UARNIHG  -  POSSIBLE  ERROR ♦  LENGTH  OF  TEXT  = 
l  '  IS  ZERO  ♦  ' 

RETURN 

C 

200  CONTINUE 
C 

ILEN  =  I 
C 

RETURN 

END 


FUNCTION  ISJZC  N1  t  N2  ) 
C 

ISIZ  *  0 
C 

DO  10  I  *  Nt  t  N2 
ISIZ  =  ISIZ  +  1 
10  CONTINUE 
C 

RETURN 

END 


' f TEXT  t 
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